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SUMMARY 


This report contains the source code and documentation for a computer program 
used to process impedance cardiography data. The cardiodynamic measures derived from 
impedance cardiography are ventricular stroke volume, cardiac output, cardiac index 
and Heather index. The program digitizes data collected from the Minnesota Impedance 
Cardiograph, electrocardiography (ECG), and respiratory cycles and then stores these 
data on hard disk. It computes the cardiodynamic functions using interactive graph- 
ics and stores the means and standard deviations of each 15-second data epoch on 
floppy disk. This software was designed on a Digital PR0380 microcomputer and used 
version 2.0 of P/OS, with (minimally) a 4-channel 16 bit analog/digital (A/D) con- 
verter. Applications software is written in Fortran 77, and uses Digital's Pro-tool 
Kit Real Time Interface Library (PRTIL), CORE Graphic Library (CGL), and laboratory 
routines. Source code can be readily modifiec to accommodate alternative detection, 
A/D conversion and interactive graphics. The object code utilizing overlays and 
multitasking has a maximum of 50 Kbytes. 


INTRODUCTION 


The Psychophysiological Research Laboratory of the Neurosciences Branch at NASA 
Ames Research Center has been engaged in a series of ground-based investigations of 
human autonomic responses to motion sickness. With the final goal of developing a 
treatment for the motion sickness-like symptoms which affect astronauts exposed to 
the microgravity environment of space, our group uses noninvasive electrophysiologi- 
cal measures to document changes in physiological activity levels in different 
subject populations. In the course of this research, we have found that measures of 
cardiovascular function are very sensitive indices of the malaise levels experienced 
by test participants. In previous studies (ref. 1), it was observed that high- 
susceptibility subjects tended to produce more labile and larger magnitude changes 
in heart rate and blood volume pulse (a relative measure of peripheral resistance), 
when exposed to motion sickness stressors than low-susceptibility subjects. We 
decided to investigate this result further. 


*University of California, San Francisco, CA. 
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Aichough electrocardiography (ECG) provides excellent information regarding the 
electrical function of the heart muscle, it gives no definitive information regard-^ 
ing physical function. Heart rate alone is insufficient for determining the level ot 
sympathetic influence on cardiac muscle. An increase in heart rate can be caused by 
either an increase in sympathetic tone or a decrease in parasympathic activation. 
Finger pulse volume, although easy to monitor, is not as reliable an index of sympa- 
thetic tone of the vasomusculature as total peripheral resistance, which is derived 
from both blood pressure and cardiac output (ref. 2). 


We needed a noninvasive, atraumatic means of examining human cardiodynamlcs in 
a moG ion- 3 icxness inducing environment. Our assessment oi impedance cardiography is 
that it provides relatively accurate data on a wide range of dynamic cardiovascular 
variables and appears even to be sensitive to rapid changes in these variables. 
Im.pedance cardiography techniques enable the measurement of myocardial contractil- 
ity, (heather index) which is directly related to the level of sympathetic innerva- 
tion of the heart. Even when there is no apparent change in heart rate, changes in 
stroke volume or myocardial contractility may be occurring which reflect significant 
differences in autoncmic function (refs. 2,3). Other valuable information obtained 
from this device include cardiac output and central volume (transthoracic 
impedance) . 


Miller and Horvath (ref. 3) describe the advantages and drawbacks of impedance 
cardiography and compare the accuracy of this technique to other invasive and nonin 
vasive measures of cardiodynamics. The principal disadvantage is that erroneous 
results may be produced if the user is unfamiliar with the effects of artifact on 
the computation of specific data epochs and is not careful to select cardiac cycles 
(dZ/dT peaks) that occur only during the expiratory pause between breaths. The 
Minnesota Impedance Cardiograph, manufactured by Instrumentation for Medicine, Inc. 
Greenwich, Conn., is available with a (firmware only) microprocessor which is 
designed to automate this process. However, we found that with the microprocessor, 
one can obtain data from one heartbeat at a minimum of 8-sec intervals. Thus much 
data are lost from the intervening beats, especially when rapid changes are occur- 
ring in the cardiovascular system. 


This paper contains the source code and operator's instructions for a user- 
interactive program written in support of our research. The program provides more 
accurate calculations of cardiac parameters based on a greater quantity of data. 

The sampling rate for digitizing data is 200 samples/sec; the data can be sampled at 
twice real-time speed from analog tape. The collected data consist of the first 
derivative of pulsatile thoracic impedance change (dZ/dt), basal impedance, electro- 
cardiogram, and respiration waveforms. During an interactive graphics session, data 
for dZ/dt, ECG and chest respiration are displayed on a monitor in 15-sec epochs. 

The user can select from the screen the dZ/dt peaks that are used to calculate 
cardiodynamic functions. The program computes and stores the means and standard 
deviations for the cardiac measures for each 15-sec epoch. 
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CALCULATIOWS PERFORMED: 


Stroke Volume 


PL'^ T(dZ/db)mLn 
. 2 


P = 53.2 

Expressed in ohms-crn 
e = 2.7133 the natural exponent 
Hematocrit expressed i i units of percent 
L = distance between t;o inner electrodes in cm 
Zo = Transthoracic Impedance expressed in ohms 


Cardiac Output = (Stroke Volume) { Pulse rate)/1000 
Expressed in liters per min 
Pulse rate is expressed in beats per min 


Cardiac Index = Cardiac Output/Body Surface Area 
Expressed in liters per min 
Body surface area = H x W x 0.00718^ 
Expressed in m^ 

H = height in cm 
W = weight in Kg 


3 


1 -eather Index - ( dZ/dt )m in/' R-Z 

2 

Expressed in cams per sec 

H-Z = inder'val between the R-wave of the EGG an 

the peak of dZ/dt expressed in sec 


REQuIRiiD ;iAhs«'jeRii AhD subJiieT Iwr’jL'.i'li-i i lOi'-j 


Hardware for this research includes a Minnesota Impedance Cardiograph, a respi- 
ration transducer (e.g., a piezoelectric or mercury strain gage) and a preampliiier 
capable of producing an analog output signal of respiratory responses, electrocar d i- 
ography equipment (eitner the Minnesota Lmpeaance cardiograph lor direct measurement 


; L I -L -rf : 


for external me?,o'.;r-=ment) , a Digital PR03S0 microcomputer wnicn 
uses version 2.0 of P/OS, and a four-channel A/D converter (l6-bit). Applications 
software, V'/ritten in Fortran 77, uses Digital's Pro-tool Kit Real Time Interi.-ce 
Library (PRTIL), CORE Graphic Library (CGL), and laboratory routines. Source cede 
can be readily modified to accommodate alternative routines for peak detection, A/D 
conversion and interactive graphics. The object code uses overlays and multitasking 
and has a maximum of 50 Kbytes. 


Additional subject information required to implement this software is 
(a) hematocrit count; (b) weight in Kg; (c) height in cm; and (d) the distance 
between the two inner impedance cardiography electrodes (tapes), measured both in 
front and back and then averaged. 


This software is installed on a DEC PRO 380 computer by following applications 
installation instructions in Professional Developer's Toolkit Refe rence Manual, 
Chapters 3 and 6. 


OftiGJNAL PAGE 

OF POOR QUALITY 
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OPERATOR IMSTRL'CTIOMS 


Three options will appear on the screen: 


* 

* SELECT FROM THE FOLLOWIMG 

* 

^ I. DIGITIZE A DATA FILE FROM TAPE 
^ 2. DATA REDUCTION OF DIGITIZED FILE 

* 3. exit 




[ENTER 1, 2, OR 3]: 

Procedure for digitizing a tape data file — SELECTION 1. If the user selects option 
1, the folicuing displays appear on the screen: 

********************* i****************** 


DIGITIZING A DATA FILE 

***********»*»***»*******»»*»***»»****»* 


PLEASE ENTER RUN ID FOR THIS FILE: 

1. User enters an ID of 6 digits. 

PLEASE ENTER RUN TIME (IN SECS.) FOR TAPE FLE: 

2. Determine maximum duration of data file in seconds (e.g. 30 minutes=l800 
seconds); ADD 60 seconds to run length (e.g , i860). ENTER THIS NUMBER. NOTE: an 
additional 60 seconds is added to the file Length to accommodate acquisition of 
calibration data. 


HOW FAST TO SAMPLE '"HE DATA ON TAPE? 


ENTER 1 FOR REAl.-TIME 


ENTER 2 FOR TWICE REAL TIME 


: ENTER 1 OR 2]: 
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3. The option of digitizin'^ data at tviloe real-time is determined by the memory 
capacity (512 K RAM) and the I/O response time (less then 0.025 sec) of the system 
implementing this program. This I/O response time is required to store four- 
channels of data^ in two-byte segments. After entering 1 or 2, the scre^^n displays; 

HIT TO START, "P" TO PAUSE, OR ”A’^ TO ABORT: 

FOLLOWED BY A [RETURN] 


4. Position analog tape to beginning of run, then press the keys "S'’ and 
"Return" to begin computer acquisition of high calibration data. Data aquisitlon may 
be paused at any time by pressing the key "P" and will not continue until the user 
presses the "RESUME" key. Acquisition may be aborted at any time by pressing the 
key "A", at which time the program prompts to the user to either SAVE or DELETE the 
created digitized file. If the user has chosen to start data acquisition (i.e., have 
pressed the key "S"), the screen will display a flashing message and three columns 
will scroll continuously until acquisition ends. The first column indicates the 


active buffer (switch*; 


between buffer 


and 


The second column in'd’" cates er^^'C'^ 


status (e.g., if A/D is turned off). Mo error is indicated by "0", 
column displays the data second currently being acquired. 


The third 


ORSGt'NA-.L , : , 

OF FCO.?? QUALH'r 
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* DATA ACQUISITION IN PROGRESS * 

1 0 1 

1 0 1 

1 0 1 

1 0 1 

1 0 1 

2 0 2 

2 0 2 

2 0 2 

* • • 

• • • 

• • , 

1 0 30 

1 0 30 

1 0 30 

IMPACQ — PAUSE - ADVANCE TAPE TO START OF LOW CALS NOW, 

THEN PRESS RESUME TO CONTINUE 

5. After 30 seconds of high calibration data, the program will automatically 
pause. Advance tape to start of low calibration data. Press the key "RESUME". 

Again, the display indicates that data acquisition is in progress. 

6. After 30 seconds of low calibration data, the program will automatically 
pause. When tape has been positioned at star:; of the data session, computer acqui- 
sition is started by pressing the key "RESUME". 

7. When the total number of seconds defining this file have elapsed the program 
screen display will direct the user to: SAVE the file; or DELETE the created digi- 
tized file. The screen displays: 
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SELECT ONE OF THE FOLLOWING [ENTER 1 OR 2] 


1 TO SAVE DIGITIZED FILE 


2 TO DELETE DIGITIZED FILE 

8. If the user has chosen to SAVE the file, it is written to hard disk, and the 
user- may now proceed with data reduction to calculate stroke volume, cardiac outpu , 
cardiac index and Heather index. The original menu is displayed: 

* SELECT FROM THE FOLLOWING 

* 1. DIGITIZE A DATA FILE FROM TAPE 

* 2. DATA REDUCTION OF DIGITIZED FILE 

» 3. EXIT 


[ENTER 1 , 2, OR 3] : 


Procedure for data reduction — SELECTION 2 

1. The screen will request input from the user as follows: 


PHASE 1 — KEY FIELD ENTRY 
a ****************************** 


ENTER RUN ID FOR THIS FILE: (6 digit file name used when digitizing). 
ENTER LENGTH (CM) BETWEEN INNER ELECTRODES: (e.g., 26.75) 

ENTER HEMATOCRIT COUNT: (e.g., 40) 

ENTER HEIGHT (CM): (e.g., 183) 

ENTER WEIGHT (KG): (e.g., 94.5) 

PHASE 2 -- CALIBRATION VALUES 
^^ssjfsaasaas************************************* 

2. The display will show the calibration levels set internal to the program 
which the user may choose to modify. These default values are: 
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Low 


High 

BASE IMPEDANCE (OHMS) 25.0000 0.0000 

dZ/dT (OHMS) 1 0000 0,0000 

WOULD YOU LIKE TO CHANGE THESE VALUES (Y/N) 

NOTE: If the user enters "N'‘, the screen will display: 

PHASE 3 -- CALIBRATION ACQUISITION 


CALIBRATION OF BASELINE IMPEDANCE 


CHANNEL 

SIGNAL 

A/DHIGH 

1 

IMPEDANCE 

5400 


DO YOU WANT TO RE-RUN HIGH CALS [Y/N]? 

NOTE: If the user enters "Y" then the the program will return to the original menu 
and user must redigitize calibration data from tape. The condition that would 
require a "Y" response here is obtaining an A/DHIGH value which does not correspond 
to the voltage out of the Minnesota Cardiograph's internal "Hi CAL" setting. In this 
example, an A/D value of 5400 equals 0.8 volts. For additional information on deter- 
mining the ratio of A/D values to voltage, refer to the: Pro/Tool Kit Real-Time 
Interface Library (PRTIL) User's Guide . Chapter 7, p. 12. 

If the user enters "N" then the screen displays 


CHANNEL 

SIGNAL. 

A/DLOW 

1 

IMPEDANCE 

433 


DO YOU WANT TO RE-RUN LOW CALS [Y/N]? 


NOTE: If the user enters "Y" then the the program will return to the original menu 
and user must redigitize calibration data from tape. If the user enters "N" then 
the screen displays: 
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S IGMAL 

ALLOW 

^DHI 

LCWCAL 

H " C A L 

1 

IMPEDANCE 

433 

5 ^^00 

0.0000 

25.0000 


MOTE: If the user enters ”Y" then the the 
user siust redigitlze calibration data 
the screen displays 


program will return to the original menu 
from tape. If the user enters ”M’’ then 






CALIBRATION OF dZ/dt SIGMhL 


CnHi'lNc.u 

SIGNAL 

AD LOW 

ADHI 

LOWCAL 

HICAL 

2 

dZ/dt 

491 

1495 

0.0000 

1 .0000 

c r .o ^ 0 

C; L- w r U 7 . 

-1 .-s « r— r- ) I 

. 7 C ' J ; C U “ O 

r ^C ^ P*^ 

“ 0 . 4 3 7 ' j 4 3 Ci 




WOULD YOU LIKE TO REPEAT CALIBRATIONS [Y/N]? 

NOTE: If the user enters "Y" then the program will return to the original menu and 
user must redigitize calibration data from tape. If the user enters "N" and this is 
the first time these data are to be analyzed, the program will display an informa- 
tive error message: -FILE DOES NOT EXIST-. It will then create an interim file and 
proceed directly to DATA REDUCTION AMD EDITING. If this is an editing session, the 
screen will display a list of all 15-second data epochs from which the user has 
selected at lease one waveform for processing. If no waveform was chosen by the user 
for a particular epoch, that epoch is "ZEROED" (i.e., the epoch number is not listed 
below and will ncc be used). A 10-min file where all epochs contain selected data is 
displayed as follows: 

FROM THE RECORD, THE FOLLOWING EPOCHS ARE NONZERO: 


1 

2 

3 

4 

5 

6 

7 

8 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


P.TE3S 

[RET 

URN] 

TO 1 

CONTIl 

NUE 
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ORIGINAL PAGE IS 
OE I»OOR QUALITY 


When the [RETURN] key is pressed, the screen will display: 


DATA REDUCTION AND EDITING 




ENTER STARTING TIME CODE FOR DATA ON TAPE (e.g., enter 300 for 3:00). 

ENTER EPOCH FOR START OF DATA REDUCTION: (e.g., 1) 

ENTER LAST DATA EPOCH: (e.g., 40) 

The time code entry is recordkeeping information which is not used by the 
program. It is here for the convenience of the user so that he may keep track of 
where in a particular tape file data reduc:;ion was started. If this is an editing 
session, the user may choose to enter only those epochs needing correction, i.e. 
new waveforems will be selected. At this point, the screen will display the data of 
the first 15-sec epoch selected: 



Figure 


1 .- 


Example of a screen display showing a 


15-sec data epoch. 
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The first waveform displayed in this figure is the dZ/dt. The horizontal lines 
through this waveform represent high cal level (upper line) and the 1°^ cal level 
(lower line). The two horizontal dots beneath each dZ/dt waveform mark the poin 
where the rising (major peak) dZ/dt waveform first crosses the low-cal line and e 
maximum trough (downward spike) following this peak. The distance between these 
dots (for a single waveform) is the ventricular ejection time. 

Th<- numerics displayed beneath the dZ/dt identify this epoch nimiber and time 
code at the start of this epoch. In this figure, the display 1 300 means epoch 1, 
time cocie 3*00. 

The second waveform is the electrocardiogram (ECG) and the horizontal dots 
beneath the ECG mark the peak of the R-wave. The third waveform shows respiratory 
cycles (Note that this waveform is inverted, such that the top of each cycle is 
played actually represents the expiratory pause between breaths). Vertical das e 
lines denote the dZ/dt peaks and their relation to ECG and respiration signa s. ac 
dZ/dt peak (vertical dashed lines) is numbered at the bottom of the screen. 

IVhile this display is on the screen, the user is prompted with questions (dis- 
played in upper left corner). The first question is. 

SKIP THIS EPOCH? [ENTER Y/N]: 

If the user does not want to have this data epoch represented in the stored 
output file, entering "Y" would cause the program to display the second 15-sec da a 
epoch. No data from the first epoch are stored, and the second epoch is renumbere 
as epoch 1. If the user enters "N" in response to this question, the next query 

displayed is 

REDO ANALYSIS? [ENTER Y/N]: 

One condition that would make the user choose to redo the analysis is that the 
dZ/dt waveform is significantly distorted by artifact (e.g., electrical noise or 
moveinent) and tharafore the horizontal dots banaaths tha wayaform ara J-orractly 
posltlonad (l.e., ventricular ejaotlon tinia Is Incorrectly labeled). Ansuering 
to this question results in the prompt 

ENTER GATE AND DIFFERENCE FACTORS SEPARATED WITH A COMMA: 

NOTE- The GATE and DIFFERENCE factors are indigenous to the peak detection subrou- 
tine used in this program. They are used for determining the peak and delta-t (ven- 
tricular ejection time) values of the dZ/dt waveform. Both factors provide a crite- 
rion for distinguishing between high frequency noise and significant trends m t e 
waveform. The GATE factor represents the minimum number of consecutive da ^ ^ 

(A/D values) needed to establish a trend. The DIFFERENCE factor represents the 
minimum difference in magnitude between data points that is needed to establish a 
real increase in the trend. If the current maximum A/D value is greater than the 
current maximum plus the minimum DIFFERENCE (current minus previous maximum A/D 
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value), then the counter is incremented. Wh:?n the counter is equal to the GATE, a 
trend toward significant peak is escabiishe'i . The default value for GATE is 2 and 
for DIFFEREMCE is 100. If the user feels that the points marking delta-t are ’’too 
close together^', then a higher DIFFERENCE factor is entered (e.g., 2,200 or 2,300). 
If the delta-t marks are ”too far apart**, t'len a lower DIFFERENCE factor is entered 
(e.g., 2,50 or 2,75). The GATE factor is rarely modified. After new GATE or DIFFER- 
ENCE factors have been entered, the screen olanks, and this epoch is again displayed 
with the dots marking delta-t repositioned. This process can be repeated until the 
user is satisfied that most dZ/dt waveforms in this epoch have delta-t correctly 
marked. If the user enuers "N**, the next qu-?ry displayed is 

OK TO CALCULATE HEATHER INDEX? [ENTER Y/N]: 

If the user enters **N*’, che output filr? would contain zeroes for the heather 
index. The condition that would make the us :*r choose to zero the heather index is 
when ^.he ECG wa\^e^orm is so distorted by artifact that the peak of the R-wave cannot 
be aetected (i.e., marked by a dot) by the ::)rogram. If the user enters '*Y**, the out- 
put file will contain the calculated heaths • inde:c fcr those waveform peaks selected 
below. The next instruction from the prograjri is 

LIST SELECTED PEAKS, SEPARATE WIT i COMMA: (e.g., 1,5,9,10,14,15) 

The criteria for selecting a peak are: (1) no significant artifact should be 

present in this waveform (and in the waveform immediately preceeding it), 

(2) delta-t must be correctly marked, and (3) waveforms must occur during the 
exhalation respiratory plateau (pause between breaths). Refering to figure 1, the 
user would likely select peaks 1,5,9,10,14, and 15. 

After making peak selections, the user presses the **RETURN** key and is prompted 
with the query 


DO YOU WANT A COPY [Y/N]: 

Answering **Y** will produce a hard copy of tne epoch being displayed on the screen. 
Answering "N** will prompt the query 

DO YOU WANT TO EXIT [Y/N]: 

Entering **N’* causes the program to display ::he next 15-sec epoch of data from which 
the user may select peaks for data reduction. This is repeated until all 15-sec 
data epochs of this file have been processed. Note that if the user has just com- 
pleted the last epoch of data in this file, the program terminates automatically 
(i.e., responds as though a ”Y” was entered). Entering ”Y’* terminates the data 
reduction phase of the program and the screen will display options for user 
selection: 
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ORIGJNAL PAGE fs 
OE ttOOR QUALITY 



EiMTER 1 TO SAVE DATA OW FLOPPY AND EXIT 


ENTER 2 TO EXIT (SAVE DATA ON HARD DISK) 

ENTER 3 TO RESTART DATA EDITNG 

ENTER 4 TO RESTART CALIBRATION ACQUISITION 

[ENTER 1, 2, 3 OR 4]: 

Entering ”1" will cause the program to prompt the user for additional key field 
information used to identify this file. The screen displays 

GnOu? TYPE (maximum of 8 characters); (e.g., CONTROL) 

RUN NUMBER: (e.g., 1) 

DIRECTION (cc,cw,nd): NOTE: cc=counterclockwise, cw=clockwise 

nd=no direction (i.e., no rotation) 
SUBJECT'S INITIALS (first and last): (e.g., PC) 

AGE: (e.g., 49) 

SUSCEPTIBILITY (1,2,3): NOTE: 1 = high motion sickness 

susceptible, 2 = moderate and 
3 = low. 

SEX (m/f): (e.g., M) 

TEST DATE (mmmddyy): (e.g., MAY0288) 

TEST TIME (military in hr, min, e.g., for 1:00 pm enter 1300): 

PRE HEART RATE (in beats/min): (e.g., 68) 

PRE TEMPERATURE: (e.g., 97.6) 

PRE B.P. (sys,dia, e.g., 120,80): 

PRE/POST TEST BASELINE (minutes): (e.g., 10) 
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0E P-OOR QUAliTY 

PLEASE IM3ERT DATA FLOPPY WITH AT LEAST 60 CONTIGUOUS FREE BLOCKS 
INTO DRIVES DZ1 & DZ2. WHEN READY, PRESS RETURN 

ARE YOU FINISHED WITH DIGITIZING )ATA? [Y/N] 

Entering "Y" will prompt the user as follows: 

DO YOU WANT TO DELETE THE DIGITIZED FILE? [Y/N]: 

(NOTE: After this input, the program returns 

the user to the original menu). 


Entering "2" will return the user to the original menu for this 
program . 

Entering "3" will return the user to the editing portion of 
this program. 

Entering "4" will return the user to the calibration acquisition 
portion of this program. 
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APPEMLIX 
SOURCE i:0DE 


PROGRAM IMPMLT 

AN IMPEDANCE CARDIOGRAPH PROGRAM USING DIGITIZED DATA FILES 
PARAMETER (IFREQ=200) 

DIMENSION SL0PE(2) ,ENTRCP(2) ,CALVAL(2,2) 

COMMON ISPACE(4^IFREQ ) , SPACE( 30*IFREQ ) 

REAL*8 DRAW ,IMPACQ , IMPCAL , IMPSTR 


DATA ACQUISITION OF BIOFEEDBACK PARAMETERS 


PEAL TEMPRT, 

i ZERO, CFREQ,TFREQ, EPOCH 

INTEGER PRED I A , PREHR , PRETEM , PRESYS , P3TD I A , PSTSYS , NUMDI A , 

@ PSTHR , PSTTEM , RUNNUM , I DA Y , NSN , AGE , SUSC , 

@ TEST I M , MI N , PSTTST , P A RT , CHANEL , TOTMIN , PRETST , 

@ IMINS,IANS, HIGH( 2,2) ,DIAG(15), 

@ CLEANS(2),ICLEAN(2) 

LOGICAL* 1 SEX , ANS , TEMP , DMY , I RSP , I ABO RT , I PAWS , TYMEAN , STDDEV , 

@ GROUP(8) ,DATE(7) ,DIRECT(2) ,03N(2) ,STRNG(8) , 

§ FRMT(12) 

CHARACTER*! CMT(80) ,RUN(4) ,RUNTYP,DATFIL( 14) ,EXP 
CHARACTER*i4 CATFIL 
EQUIVALENCE (DATFIL, CATFIL) 


INITIALIZATION 


DATA IMPCAL/ 'IMPCAL'/ , IMPSTR/ ' IMPSTR' / 

DATA DRAW/'DRAWRM'/ , ICLEAN/27,99/ , IMPACQ/ ' IMPACQ' / 

DATA CLEANS/27,99 /,CATFIL/'00000000000000'/ 

TYPE 500, CLEANS 
500 F0RMAT(X,4A1) 

TYPE 101 

101 F0RMAT(10X,60('*' ) ,//,27X, ' IMPEDANCE CARDIOGRAPH PROGRAM' , 

1 //,10X,60( '*') ,/////) 

125 CONTINUE 

CATFIL( 1:14)=' 00000000000000 ' 

TYPE 102 
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' 02 F CRMAT( 20X . 45 . 20.’( , ' * ' , 

1 / , 20X 6X . ' SELECT FROM THE FOLLOWING ’ , / , 20X , ' * ' , 

1 /,20X, , 10X, ' 1 . DIGITIZE A DATA FILE FROM TAPE’, 

2 /,20X, , 10X, '2. DATA REDUCTION OF DIGITIZED FILE’, 

3 /,20X, , lOX, ’3. EXIT' . 

4 /,20X,45(’-’),/, 

5 /' [ENTER 1, 2 OR 3J: '$) 

150 CONTINUE 

ACCEPT * ,IAMS 
IFdANS.EQ. 1) GOTO 200 
IF(IANS.EQ.2) GOTO 2000 
IF( IANS.EQ.3) GOTO 99999 

TYPE *,’ *»*^*^+ ERROR- TYPE 1 OR 2 OR 3 ONLY. NOW DO IT' 
GOTO 150 
200 CONTINUE 

IEFN=6 

C CALL IMPACQ(CATFIL) 

CALL SPAV/N{ RAD50( IMP ACQ ♦ lEFM , I ESD ,,,,,, IDS ) 

D TYPE SPAWN CALLED ' , ' IEFN= ' JEFN, ' IESD= ’,IESD,'IDS= ’ , IDS 

CALL WAITFR(IEFN,IDS2) 

C TYPE lEFN SET NOW ’,IDS2 

TYPE 500, CLEANS 
GOTO 125 
2000 CONTINUE 

C ========================== 

C PHASE 1 — KEY FIELD ENTRY 

C ========================== 

TYPE 500, CLEANS 
TYPE 501 

501 FORMAT( 10X,60( ’*'),/////, 27X, ' PHASE 1 — KEY FIELD ENTRY',/////, 
1 10X,60( ’*’) ,/////) 

TYPE 3998 

3998 FORMAT (/' ENTER RUN ID FOR THIS FILE:’,$) 

ACCEPT 3999,CATFIL(5: 10) 

3999 F0RMAT(a6) 

TYPE 4000 

4000 FORMAT(/’ LENGTH (CM.) BETWEEN INNER ELECTRODES: ',$) 

ACCEPT » ,ELEN 

TYPE 4100 

4100 FORMAT(/' HERMATOCRIT COUNT : ',$) 

ACCEPT * ,HEM 
TYPE 4200 

4200 FORMAT(/’ HEIGHT(CM) : ',$) 

ACCEPT * , HEIGHT 

TYPE 4300 

4300 FORMAT(/' WEIGHT(KG) : ',$) 

ACCEPT » , WEIGHT 


SET DEFAULTS AND ZERO CURRENTLY UNUSED FIELDS 

original page fS 

OF POOR QUALITY 
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RHO=53.2*EXP (.022*HEM) 

BSA=HE IGHT^» . ?25^WE IGHT*» , 425* . 007 U!4 
C=RHO*ELEN**2 

D TYPE ' C= ' ,C, 'RH0= ' , RHO 

D CALL WAIT(5,2) 


C PHASE 2 — CALIBRATION VALUES 

6000 CONTINUE 

CALL ERRSET{59, .TRUE. , .FALSE. , .TRUE. . .FALSE. .MAX) 
CALL IMPVAL (CALVAL) 

C PHASE 3 -- CALIBRATION ACQUISITION 


7010 CONTINUE 

OPEN (UMIT=1 ,MAME: 'DU1 :NUDRAW.TEL' ,T:'PE='MEW , FORM= ' FORMATTED ' ) 
WRITE( 1 ,2) CALVAL ,C, BSA,CATFIL , IMINS 

2 F0RMAT(6( G15.7),A14,I10) 

CLOSE (UNIT=1) 

IEFN=6 

CALL SPAWN (RAD50(IMPCAL),,,IEFN,,It:SD,,,,,,IDS) 

D TYPE SPAWN CALLED ' , ’ IEFN= ',IEFN„' IESD= ',IESD ,'IDS= '.IDS 

CALL WAITFR(IEFN,IDS2) 

C TYPE lEFN SET NOW ',IDS2 

OPEN (UNIT= 1 ,NAME= ' DW1 rNUDRAW.TEL ' ,TYPE= 'OLD' , FORM= ' FORMATTED ' ) 
READ (1,2) SLOPE, ENTRCP,C,BSA,CATFIL, IMINS 
CLOSE (UNIT=1) 

CATFIL(11:14)=' .IMP' 

OPEN INTERNAL FILE THAT STORES RESULTS OF THE REDUCTION PROCESS 
IF IT EXISTS AND FIND OUT WHAT EPOCHS HAVE BEEN COMPLETED. OTHER 
WISE CREATE A NEW INTERNAL FILE TO BEGIN REDUCTION 

OPEN (UNIT=1,STATUS='0LD' ,NAME=CATFIL, 

1 F0RM=' UNFORMATTED' , ACCESS= ' DIRECT' , REGL=8, ERR=50) 

CL0SE(UNIT=1 ) 

CALL INVEST(CATFIL) 'DETERMINE WHAT EPOCHS HAVE BEEN COMPLETED 
TYPE *, 'PRESS [RETURN] TO CONTINUE ' 

ACCEPT 500 ,ANS 
GOTO 55 
50 CONTINUE 

OPEN (UNIT=1 ,STATUS='NEW' ,NAME=CATFIL , 

1 FORM=' UNFORMATTED' , ACCESS= ' DIRECT' , RECL=8) 

CL0SE(UNIT=1) 

MAXREC=380 

IF(CATFIL( 10: 10) .EQ. 'O' ) MAXREC=120 
IF(CATFIL( 10:10) .EQ. '4' ) MAXREC=220 

CALL ZEREC (MAXREC,CATFIL) IZERO CONTENTS OF NEW FILE FOR MAX RECS 
55 CONTINUE 


19 



CO o 



9101 CONTINUE 
TYPE 9109 

9109 F0RMAT( //////////) 

TYPE 9110 

9110 FORMAT ( ' WOULD YOU CARE TO 

1 /7X/1) SAVE DATA ON FLOPPY AMD EXIT', 

2 /7X,'2) EXIT (SAVE DATA OM DISK)', 

3 /7X,’3) RESTART EDITING', 

4 /7X,'4) RESTART CALIBRATION ACQUISITION (SAVE)', 

5 /' [ENTER 1 ,2, 3, OR 4 ] : ' ,$) 

ACCEPT 10, TEMP 

10 FORMAT(AI) 

IF(TEMP.EQ. '2' .OR .TEMP . EQ . ' 4 ' ) THEN 
OPEN{UNIT= 1 , NAME= 'DW1 :NUDRAW.TEL' ,STATUS= 'OLD' ) 

CLOSE (UMIT=1,DISP0SE='DELETE' ) 

END IF „ .X 

IF (TEMP .EQ . ' 1 ' .OR .TEMP . EQ. ' 2 ' .OR.TEMP .EQ. ' 3' -OR .TEMP .EQ. 4 ) 

1 GO TO 9114 

TYPE 9 1 08 

9108 FORMAT(////' YOU MUST NOT HAVE CHOSEN A VALID ALTERNATIVE;'/, 

@ ' CHOOSE ONE OF THE FOLLOWING; '/) 

GOTO 9101 
9114 CONTINUE 

DECODE (1,1 10, TEMP) ITEMP 
110 FORMAT(II) 

GOTO (91 13,95990, 30000,6000) ITEMP 
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c 

C ASSEMBLE DATAFILE WAME 

C 


OE fipOR QUALITY 


9115 COMTIMUE 

CALL SPAWN(RAD50( IMPSTR) , , , lEFM, , lESD, , , , , , IDS) 

D TYPE SPA'^ CALLED ' , ' IEFN= ',IEFN,' IESD= ',IESD,’IDS= *,IDS 

CALL WAITFRdEFM, IDS2) 

C TYPE lEFM SET NOW ',IDS2 

C CALL IMPSTR( CATFIL , IFLAG , GROUP , RUNNUM,DIRECT,OSN , AGE ,SUSC , SEX , 

C 1 DATE , TESTIM , PREHR , TEMPRT , PRETEM , PRES YS , PREDI A , PRETST , PSTTST ) 

C TYPE 500, CLEANS 

0PEN(UNIT=1 , NAMEr'DWI :NUDRAW.TEL' ,STATUS='OLD' ) 

CLOSE (UNIT=1,DISP0SE= 'DELETE') 

GOTO 125 

Cl 0000 CONTINUE 


c 

IFLAG: 1 ! WRITE COMMENTARY FILE 

c 

CALL COMEMT (KCHAR) 

c 

GOTO 9101 

99990 

CONTINUE 


TYPE 500, CLEANS 


GOTO 125 

99999 

CONTINUE 


END 

SUBROUTINE ZEREC( MAXREC , CATFIL) 
integer tcount 
CHARACTER* 14 CATFIL 

OPEN (UMIT=1 , STATUS: 'OLD' ,NAME=CATFIL , 

1 FORM: 'UNFORMATTED' , ACCESS: ' DIRECT' ,RECL:8) 

do 100 tcount : 1, MAXREC 

write( 1 ,REC:TC0UNT,ERR:999)0. ,0. ,0. ,0. ,0. ,0. ,0. ,0. 

100 continue 

999 CL0SE(UNIT:1, STATUS: 'SAVE' ) 

RETURN 

END 

SUBROUTINE INVEST( CATFIL) 

DIMENSION BUF{8,1),N0NZER(400) 

CHARACTER* 14 CATFIL 

CALL ERRSET ( 36 , . TRUE . , . FALSE . , . TRUE . , . FALSE . ) 

IREC=1 

INUM:0 

TYPE *,' FROM THE RECORD, THE FOLLOWING EPOCHS ARE NONZERO :' 
100 CONTINUE 

OPEN (UNIT:1, status: 'old' , NAME:CATFIL,ERR:999 , 

1 FORM: ' UNFORMATTED ' , access='direct' ,recl:8 ) 

read( 1 ,REC:IREC ,ERR=999)(BU-( ICOUNT, 1 ) , ICOUMT:! ,8) 
IF( BUF(1,D) 999,999,200 

200 CONTINUE 

INUM=INUM+1 

NONZER(INUM):IREC 

999 CL0SE(UNIT:1 , STATUS: ' SAVE ' ) 
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'.’ALL ZLAL'L T ' Lo , I ZRR ) 

IF(IERR.EQ. ! ) THEM 

TYPE *,(M0MZER(I),I=1,IMUM} 

Q r^iTir f p \ j 

ELSE 

IREC^iRECt 1 
GOTO 100 
END IF 
EMD 

PROGRAM IMPCAL 
C (SLOPE, IMTRCP,CATFIL) 

PARAMETER ( IFSEQ^EOO; 

REAL CALVAL{2, 2) , INTRCP(2) ,SL0PE{2) 


C 

C 


2 

7010 

500 

700 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


c 

c 


INTEGER CLEAMS(2) , CHANEL, HIGH{2, 2) 
COMMON IS?ACE{4^IFREQ) ,SPACE( SO^IFREQj 
EQUIVALEMCE(CALVAL,SPACE( 1 ) ) 

''''u i P !i r’A'T'rTf 

DATA CLEANS/ 27, 99/ 


PHASE 3 — CALIBRATION ACQUISITION 


Ci* RiiQ= If REQ _ 

OPEN(UNIT= 1 ,NAME= ' DW i iMUDRAW.TEL' ,TYPE= 'OLD' ,ORM= ' FORMAiTbD ' ) 
READ (1,2)CALVAL ,C,BSA,CATFIL, IMINS 

F0RMAT(6(G15.7),A14,I10) 

CL0SE(UNIT=1 ) 

TYPE 500, CLEANS 
F0RMAT(X,4A1) 

TYPE 700 

FORMAT(10X,60( '*'),/////, 22X, ' PHASE 3 — CALIBRATION 
ACQUISITION' ,/////, 1 0X, 60( ' ),/////) 

CATFIL( 1 ;4)-'DWl : ' 

CATFIL( 11:14)=' .DIG' 


CALIBRATE BASELINE IMPEDANCE - CHANNEL 2 


CALL CALZ (SL0PE(1),INTRCP(1),CALVAL,CATFIL) 


CALIBRATE dZ/dT SIGNAL - CHANNEL 3 


CALL CALDZT (SL0PE(2) , INTRCP(2) ,CALVAL,CATFIL) 

0PEN(UNIT=1 ,NAME='DW1 iNUDRAW.TEL' ,TYPE='OLD' , FORM= ' FORMATTED ’ ) 
WRITE( 1 ,2)SL0PE,INTRCP,C,BSA,CATFIL, IMINS 
CL0SE(UNIT=1 ) 

END 

SUBROUTINE DIRECT (TCOUNT,CATFIL,NCHAN) 

STORAGE SEQUENCE AS FOLLOWS: 

FIRST - RAW RESPIRATION ^ 

OF POOR QUALITY 
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C SECOMD - DELTA Z QUALITY 

C THIRD - DZ/DT 

C FOURTH - ECG 

c TKB FOR TASK BUILD 

u 

c TKB>FILENAME=FILEMAME 

C TKB>/ 

C TKB>MAXBUF=1600 

c TKB>// 


TOO 


PARAMETER (IFREQ=200) 
integer buf( IFREQ ), tcount 
DIMENSION X( 15*IFREQ) , Y( 15*IFREQ) 
COMMON BUF,X,Y 
CHARACTER*14 CATFIL 


OPEN (UNIT=1, status='old' , NAME=CATFIL 

•D' 


,recl=4*IFREQ) 


read( K REC=TCOUNT, ERR=999) (BUF( ICCUNT ) , ICOUNT= 1 , 4»IFREQ) 
continue 


999 CLOSE(UNIT=1 , STATUS='SAVE' ) 

10 FORMAT(x,i3,4(2X,l4,2x,i4)) 

20 FORMAT( 1X,I5) 

RETURN 
END 

SUBROUTINE CALZ (SLOPE, INTRCP,CALVAL, CATFIL) 
CALIBRATE BASELINE IMPEDANCE - CHANNE,. 2 


REAL CALVAL{2, 2) ,CFREQ,TFREQ,EPOCH, INTRCP 

INTEGER*4 lADSUM 

INTEGER CHANEL, ICOUNT,HIGH( 2), 

@ CLEANS ( 2 ) , MAXPTS , 1ST AT ( 2 ) , lEFN , MDSYM , 

@ ICHAN(8) , ICONV, IFORM, ITRIG, ISTAT2(2) ,IEFN2,NPTS(2) 

CHARACTER»4,IWORD(2) 

CHARACTER»14, CATFIL 

PARAMETER (IFREQ=200) 

COMMON ISPACE{4*IFREQ),SPACE(30*IFREQ) 

LOGICAL* 1 ANS 


DATA IWORD /'HIGH', 'LOW '/ 
DATA CLEANS/27,99/ 

C 

C 

c 
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INITIALIZATIGM 


0 

C 

c 

c 

c 

c 


500 

TOO 


PERFORM PRELIMINARY 


ACQUISITION 


CALL WAIT( 2,2) 

TYPE 500, CLEANS 
F0RMAT(X,4A1 ) 

TYPE 700 

F0RMAT(10X,60( ' ' ),/////. 22X 
OF BASELINE IMPEDANCE ' . ///// 


' CALIBRATION 
10X,60( )//////) 


CALL WAIT( 2.2) 

CHANEL=2 

C 

C 

7001 CONTINUE 

r 

C 

C 

c 

C 0PEN(UNIT=2 ,NAME= ' SY :DUMP .TST ' , FOftM= ' FORMATTED' ,TYP£= 'NEW ) 

ISTREC=1 
7019 CONTINUE 

N=1 !HIGH CALIBRATIONS 
C 

C PERFORM HIGH CALIBRATION 

C 

C CALIB=.TRUE. 

7110 CONTINUE 

C 

C INITIALIZE SLOPE, INTERCEPT 

C 

SLOPE =0.0 
INTRCP =0.0 
NSAMP=0 
IADSUM=0 

DO 7000 1=0,29 

IREC= 30*(N-1 )+ISTREC+I 

CALL DIRECT(IREC,CATFIL) 

C 

c 

c 

c 

C CALCULATE MEAN 

C 

C 
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:i?tc:l'j=t?r£q .'number of samples fe;: ch.amnel per record=per second 

DO 704 J=1,NPTCHN 

IADSUM=IADSUM+ISPACE ( ( J- 1 ) »4+CHi,NEL ) 

0 TYPE »,IADSUM,ISPACE ( ( J- 1 ) *4^-CHt.MEL ) 

MSAMP-NSAMP+1 

C WRITE(2, 19999 ) I AD(I+(J-1 )*CHANEL+(K-1 )»MAXPTS),NPTCHN, CHANEL, 

C 1 NPTS(K),NSAMP,I,J,K,N 
C19999 F0RMAT(9I8) 

704 CONTINUE 

C WRITE(2,*) IADSUM,NSAMP,ISPACE 

7000 CONTINUE 

H IGH( M ) - 1 ADSUM/MS AMP 
C 

c 

C ECHO MEAN AND VERIFY 

C 

C 

TYPE 715,IWORD(N) 

"7 ' 5 FORMAT( , ' ■'// '//// ' CHANNEL 

1 -,5x,' 

TYPE 71550 , HIGH(N) 

71550 FORMAT( 5X,' 1 ' ,8X, ' IMPEDANCE ',5X,I10) 

TYPE 717,IWORD(N) 

717 FORMAT(/' DO YOU WANT TO RE-RUM ',A4,' CALS [Y/N]? ',$) 

ACCEPT 617,ANS 

617 FORMAT(AI) 

IF (ANS .EQ. 'Y' .OR.ANS.EQ. 'y' ) THEN 
TYPE WHAT STARTING RECORD ' 

ACCEPT *,ISTREC 
GOTO 7110 
END IF 

IF (N.EQ.2) GOTO 724 
N=2 

GOTO 7110 
724 CONTINUE 

C 

C CALCULATE SLOPE AND INTERCEPT 

C 

IF (HIGH(I) .NE. HIGH{2)) GOTO 71o5 
INTRCP =0.0 
SLOPE =0.0 
GOTO 718 

7185 SLOPE = (CALVAL(2, 1 )-CALVAL( 1 , 1 ))/FLOA': (HIGH{ 1 )-HIGH(2) ) 

INTRCP = CALVAL(1,1) - SL0PE*FL0A'''(HIGH(2) ) 

718 CONTINUE 

C 

C ECHO SLOPE, INTERCEPT 

C — 
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,5X' SIGNAL 
• 5X,- 
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,5X,' ADLOW 


TYPE 720 

720 FORMAT(/’ CHANNEL ',5X,' SIGNAL 


\5X, 


1' ADHI ',5X, 

2 ' LOWCAL ',5X,' HICAL ',5X)) 

TYPE 81550 , HIGH(2) 

31550 FORMAK 4X,’ 1 7X IMPEDANCE ',4X,I10,$) 

TYPE 722, HIGH ( 1 ) ,CALVAL(1 , 1 ) ,CALVAL (2,1) 

721 CONTINUE 

722 FORMAT( ’ + ' ,4X , 1 10 , 4X, F 10. 4 , 4X ,F 10.4) 

TYPE 723 

723 FORMAT(/' WOULD YOU LIKE TO REPEAT CALIBRATIONS [Y/N]?’,$) 
ACCEPT 617,ANS 

IF (AMS .EQ. ' Y' .OR.AMS.EQ. 'y ' ) THEN 
TYPE WHAT STARING RECORD ' 

ACCEPT »,ISTREC 
GOTO 7019 
END IF 

CLOSE (UMIT=2) 

RETURN 

END 

FUNCTION MAXIAD ( IAD,MAXPTS,MPTS) 

DIMENSION IAD(NPTS*MAXPTS) 

MAXIAD= IAD(3) 

DO 100 J=1,MAXPTS 

IF(IAD((J-1)*NPTS+3 ).LT. MAXIAD) GOTO 100 
MAXIADr IAD((J-1)*NPTS+3) 

100 CONTINUE 
RETURN 
END 

FUNCTION MINI AD ( IAD,MAXPTS ,MPTS) 

DIMENSION IAD(NPTS*MAXPTS) 

MINIAD= IAD(3) 


DO 100 J=1,MAXPTS 

IF(IAD((J-l)»NPTS+3 ).GT.MIMIAD) GOTO 100 
MINI AD= I AD( ( J- 1 ) *NPTS+3 ) 

100 CONTINUE 
RETURN 
END 

SUBROUTINE CALDZT (SLOPE , INTRCP ,CALVAL,CATFIL) 


CALIBRATE dZ/dT SIGNAL - CHANNEL 3 


REAL CALVAL(2, 2) ,CFREQ,TFREQ, EPOCH, INTRCP 
INTEGER*4 lADSUM, IAD4 , IAD4MX , I ADELT 
INTEGER CHANEL, IC0UNT,HIGH(2) , IAD(400) , 

§ CLEANS ( 2 ) , MAXPTS , 1ST AT ( 2 ) , lEFN , MDSYN , 

§ I CHAM ( 8 ) , I COMV , I FORM , ITR IG , ISTAT2 ( 2 ) , IEFN2 , NPTS ( 2 ) 

PARAMETER (IFREQ=200) 

COMMON ISPACE(4*IFREQ) ,SPACE( 30*IFREQ) 

EQUIVALENCE (lAD.ISPACE) 
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CHA.FACTER^4,IW0RD(2) 
CHARACTERS 14, CATFIL 


C 

C 


C 

C INITIALIZATION 

C 

DATA CLEANS/27,99/ 

C PERFORM PRELIMINARY ACQUISITION 

TYPE 500, CLEANS 
500 FORMAT ( X, 2A 1 ) 

TYPE 700 

700 FORMATC 10X , 60( ' s ' ) , ///// , 22X , ' CALI3RAT10M 

1 OF dZ/dT SIGNAL' ,/////, 10X,60( ' S' ) ,/////) 
CALL WAIT(5,2) 

CHANEL =3 

Q 

C 


7001 

C 

C 

C 

C 

7015 

C 

c 

c 

c 

7019 

C 


CONTINUE 


CONTINUE 


CONTINUE 

OPEN( UNIT=2 , N AME= ' SY ; DUMP . TST ' , FORMr: 


FORMATTED' ,TYPE='NEW ) 


N=1 !HIGH CALIBRATIONS 
SLOPE =0.0 
INTRCP=0.0 
713 CONTINUE 
NCHAN=0 
NSAMP=0 
IADSUM=0 


C 

c READ IN A-D VALUES 


IREC= 31+1 

TYPE IREC= ',IREC,'N= ',N,'I= ',I 
CALL DIRECT(IREC,CATFIL) 

TYPE *, ( IAD{ ( J-1 )* 4 +3),J=1,100) 
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c 

C CALCULATE MEAN 

C 

c 

I?(M.EQ. 1 ) THEM 
C EIMD MAX VALUE IF M=1 

IADMAX-MAXIAD( lAD, IFREQ ,4 ) 

D TYPE TADMAX= ' , lADMAX 

IADELT=IADMAX 
IADELT-IADELT^5 

IF( IADELT.lt. 0) IADELT=-IADELT 

C FOR M=2 FIND MAX AMD MIM 

c IADMAX-MUL( lAD.MAXPTS.MPTS) 

I .ADMAX^MAX I AD ( I AD , I F REQ , 4 ) 

TADMIM=MIMIAD{ lAD, !FREQ,4 ) 

rZEHO= I ADMIM+( I ADM AX- 1 ADMIM ) »2/5 
END IF 

NPTCHN= IFREQ 
DO 704 J=1,NPTCHM 
IF(iM.EQ.I) THEM 
C 

c USE lAD VALUE IF lAD >= MAX VALUE- 5%*(MAX VALUE) 

C 

IAD4=IAD({J-1)»4+CHAMEL ) 

IAD4MX=IADMAX 

IF(( 100^IAD4) .GE. ( IAD4MX*100-IADELT )) THEM 
MCHAM=MCHAN+1 

IADSUM=IADSUM+IAD ( ( J-1 )»4+CHAMEL ) 

C TYPE »,IAD4,IADELT, IAD4MX 

EMD IF 
ELSE 

IF(IAD((J-1)*4+CHANEL ).LE.IZERO +80.AND. 

1 IAD((J-1)*4+CHAMEL ).GE.IZERO -80) THEN 
NCHAM=NCHAM+1 

IADSUM=IADSUM+IAD ((J -1)*4+CHANEL ) 

EMD IF 
EMD IF 

C WRITE{2, 19999)IAD(J+(K-1)»MAXPTS) ,NPTCHM,CHAMEL, 

C 1 NPTS(K),NCHAM,I,J,K,IADMAX 
19999 F0RMAT(9I8) 

704 CONTINUE 

7000 CONTINUE 

HIGH{ M ) = I ADSUM/MCHAM 
IF (M.EQ.2) GOTO 714 
N=2 

GOTO 713 

714 CONTINUE 


ORIGINAL PAGE IS 
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ECHO MEANS AND VERIFY 


C 
C 
C 
C 

C WRITE(2,* ) IAD,HIGH,IADMAX,IADMIN,IZERO 

C TYPE 715,IWORD(M) 

715 FORMAT(/////////' CHANNEL ' ,5X' SIGNAL ' , 5X , ' A/D ' , A4 , 

1 /,- ') 

C TYPE 71551 , HIGH(N) 

71551 FORMAT! 5X,' 2',8X,' dZ/dT ’,5X,I10) 

C TYPE 717 

717 FORMAT!/’ DO YOU WANT TO RE-RUN dZ/dT CALS [Y/N]? ',$) 

C ACCEPT 617,ANS 

617 FORMAT! Al) 

C IF (ANS .EQ. 'Y' .OR. AMS.EQ. 'y' ) GOTO 7110 

CLOSE! UNIT=2) 

C 

C CALCULATE SLOPE AND INTERCEPT 

C 

IF (HIGH!D .NE. HIGH!2)) GOTO 7135 
INTRCP =0.0 
SLOPE =0.0 
GOTO 718 

7185 SLOPE = (CALVAL(2,2)-CALVAL! 1 ,2))/FL0AT (HIGH(1 )-HIGH(2) ) 
INTRCP = CALVAL(1,2) - SL0PE*FL0AT(HIGH(2) ) 

718 CONTINUE 

C 

C ECHO SLOPE, INTERCEPT 

C 


TYPE 720 

720 FORMAT!/' CHANNEL ',5X,' SIGNAL ',5X,' ADLOW ',5X, 

1' ADHI ',5X, 

2 ' LOWCAL ',5X,' HICAL ',/,X,6(' ’,5X)) 

TYPE 81551 , HIGH! 2) 

81551 FORMAT! 4X,' 2',7X,' dZ/dT ',4X.I10,$) 

TYPE 722, HIGH ! 1 ) , CALVAL! 1 ,2) ,CALVAL !2,2) 

722 FORMAT! ’+’ ,4X , 110 , 4X , F10 . 4 , 4X, F10.4) 

TYPE *, ' SLOPE ’, SLOPE,' INTERCEPT INTRCP 
TYPE 723 

723 FORMAT!/’ WOULD YOU LIKE TO REPEAT CALIBRATIONS [Y/N]?',$) 
ACCEPT 6 17, ANS 

IF !ANS .EQ. 'Y' .OR. AMS.EQ. 'y' ) THEN 
TYPE »,' WHAT STARING RECORD ' 

ACCEPT *,ISTREC 
GOTO 7019 
END IF 
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c 


. i ■ i J 

EMD 

PROGRAM IMPACQ 

REAL TEM?RT.CAL'/AL(2 . 2),MEAiM(2 ).SD( 2),SLOPE( 2),IMTRCP( 
3 ADSUM,VARIAM( 2 ) . 2ERO,CFREQ,TFREQ, DEV( 2), EPOCH 


2 ), 




§ 


INTEGER IOSB(2),BLIMK(3) ,RESET(2) , 
ISTA(4) ,ICONT(2) , 


MIN , CHANEL , ICLK , I ADCSR , 
ICOUMT,IMINS,IANS,LOW( 10) ,HIGH( 10) 
CLEAMS(2) ,MAXRUN,MAXPT3,ISTAT(2) ,IEFN( 
ICHaN(8) , ICOMV, IFORM, ITRIG, I3TAT2^2) , IE 
INTEGER lEFM i , IPAR{6) , I3TAT3(2) 


,DIAG( 15), 
2) ,MB3YM, 
FN2.NPTS(2) 


L 

C 

C 

0 

c 

c 


L0GICAL*1 EXP, 3EX,AN3, TEMP, DMY, TYMEAN,STDDEV , 

i? DATFIL( 14) ,STRNG(8) , 

@ FRMT(12) 

(^UApii pTcrp^lli 

CHARACTER^ 1 I8UF , IPAW3, lABORT, IRSP 

pnRAM£TEH (I^REO =200 ) ! 3ETS PROGRAM 3IZE HRT SAMPLE RATE 
PARAMETER (MDSYN=1) 

PARAMETER (MODE =0) 

PARAMETER (IC0NV=0) 

PARAMETER (IF0RM=0) 

PARAMETER (ITRIG=0) 

PARAMETER (ITIME=0) 

PARAMETER (CHANEL=4) 

INTEGER*4 lADSUM ,AD( 8),ISIZE 

INTEGER»2 IAD(8^IFREQ ) , IAD2(4*IFREQ) , IAD1 (4*IFREQ) 

EQUIVALENCE ( IAD1 , lAD) , ( I AD2 , IAD(4»IFREQ+1 ) ) , ( ISTAT, ISTA) , 

1 {ISTAT2,ISTA(3)),(IEFN(1),IEFN1),(IEFN(2), IEFN2) 

COMMON I AD 


INITIALIZATION 


500 

510 

617 

618 
7125 
70 ’7 
7135 
10 


DATA CLEANS/27,99/ 

DATA BLINK/155,53, 109/, RESET/155, 109/ 

DATA FRMT/ ' ( ’ , ' ','(','X',’,','F', '9', '.','4', ')',')'/ 

DATA ICLK/" 172540/ 

DATA IREAD/"001000/ 

DATA lABORT/ "101 /, IPAWS/"120/ 

F0RMAT(X,4A1 ) 

F0RMAT(X,5A1 ) 

FORMAT(AI) 

F0RMAT(A6) 

FORMAT(/' AIN STATUS WAS:') 

FORMATC STAT(l) (OCTAL) ^',C6,' ISTAT(2) (DECIMAL) = ’,I6) 
FORMAT ( ' TRY AGAIN. ’ ) 

F0RMAT(I2) 


ORIGSN^L Pi YI : . 
OF POOR QUALITY 
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C PHASE 5 — DATA ACQUISITION 

C =========z===z=======z===== 

FREQ=IFREQ 

MAX?TS=4»IFREO 

CALL GETADR(IPAR(1),IBUF) 

IPAR(2)=1 

9000 TYPE 500, CLEANS 

TYPE 900 

900 FORMAT( 10X,60( '*'),/////, 22X, ' DIGITIZING A DATA FILE 

1 /////, 10X, 60( '»' ),////) 

c 

C REQUEST MAXIMUM RUM DURATION 

C 

TYPE 90109 

90109 FORMAT(//,'’ PLEASE ENTER RUM ID FOR THIS DATA FILE: ',$) 

ACCEPT 618,(CATFIL(5:10)) 

CATFIL( 1 :4)z'DW1 : ' 

CATFIL( 11 :14)=' .DIG' 

TYPE 90110 

90110 FORMAT(///’ PLEASE ENTER RUM TIME (IN SECS.)FOR TAPE FILE: ’,$) 
ACCEPT *,MAXRUN 

ISIZE=6.25 *MAXRUN+1 
TYPE 90111 

90111 FORMAT(///’ HOW FAST TO SAMPLE THE DATA ON TAPE ? 

1 /10X,' ENTER 1 FOR REAL TIME RATE 

2 /I OX,' ENTER 2 FOR TWICE REAL TIME ’, 

3 /' [ENTER 1 OR 2] NOW: ' ,$) 

ACCEPT »,ISAMP 

C MOTE: MAXPTS^FREQ*4 

CFREQ=FREQ*ISAMP 
CALL SETFRQ(CFREQ,TFREQ) 

DO 190 1=1,4 
ICHAM(I)=1 
190 CONTINUE 

IEFN1=8 
IEFN2=10 
IC0UNT=2*MAXPTS 

C 

C OPEN TEMPORARY DATA FILE 

C 

9015 CONTINUE 

C 0PEN(UNIT=2,MAME='SY:TEMP.DAT' ,TYPE='NEW' ) 

OPEN (UNIT=1 ,STATUS='NEW' ,MAME=CATFIL. INITIALSIZE=-ISIZE, 

1 FORM=’ UNFORMATTED' , ACCESS= ' DIRECT RECL=4*IFREQ) 

MIN = 0 
IMINS = 0 
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C INITIALIZE 3UM,IC0UMT 

C 


9011 CONTINUE 

ICONT( 1)=IC0UNT/2 
IC0NT(2)=IC0UNT/2 
CALL AINIT (ISTAT.IEFNI) 

IF( ISTAT(1).EQ. "40000. OR. ISTAT(1).EQ. ”40040) GOTO 7000 

TYPE ISTAT(I) = ',ISTAT(D 

PAUSE 

7000 CONTINUE 


7100 

1 


9022 

9012 

@ 


97100 

1 

2 

3 


TYPE ' ' 

IF (IRSP .NE. IPAWS)GO TO 9022 
PAUSE 

TYPE 500, BLINK 
TYPE 7100 

FORMAT( / ' ir jf iTTTTT^T^TTt ^ 

/' * DATA COLLECTION RESUMED 

TYPE 500, RESET 
GO TO 901 
TYPE 9012 

FORMAT (' HIT "S" TO START, "P" TO PAUSE,’, 

' "A" TO ABORT:'/' FOLLOWED BY A [RETURN]') 

ACCEPT 617, ANS 

IF (ANS .EQ. ' A’ ) GO TO 9100 

IF (ANS .EQ. 'P') PAUSE 

TYPE 500, CLEANS 

TYPE 500, BLINK 

TYPE 97100 

FORMAT( ///////////////, 

/I ^ 

/’ » DATA COLLECTION UNDERWAY »' , 


TYPE 500, RESET 
901 CONTINUE 

CALL QI0(IREAD,7, 4 , , lOSB , IPAR , IDSW) 

CALL AIN (ISTAT, ! START CLOCK 

2 lAD , 

3 ICONT( 1 ) , 

4 IEFN1, 

5 MDSYN, 

6 ICHAN, 

7 ICONV, 

8 I FORM, 

9 ITRIG.ITIME,CFREQ) 

CALL ABUF(ISTAT2,IAD2,IC0MT(2),IEFN2,MDDE ,ITIME) 
908 CONTINUE 

IMINS = IMINS + 1 
C TYPE *, ’ IMINS= ' , IMINS 

K=1 
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o o o o n o 


C M=1 

C lENDrlCOUMT 

80001 CONTINUE 

IISTAT=ISTA(1+2'^(K-1)) 

IF (IISTAT .EQ."40000.0R.II5TAT .EQ."40040) 

1 GOTO 81300 

C IF (IISTAT. ME. 0) TYPE *,K, IISTAT 

TYPE »,K, IISTAT, IMIMS 
GOTO 80001 
81300 CONTINUE 
C 
C 

C 

c STORE DIGITIZED DATA FOR EACH SECONI' 

C 

C 

CALL DIRECT( IMIMS,K) 

NCONT=ICONT(K) 

CALL ABUF( ISTA( 1+2^(K-1 ) ) , IAD( 1+(K-1 )^MAXPTS) ,NCONT ,IEFN(K), 

1 MODE.ITIME) 


ECHO MEANS AMD STANDARD DEVIATIONS 


IF END OF SESSION, END DATA ACQUISITION 


MIN = IMINS 
IF(MIN.EQ. 30) THEM 
CALL ACSTAT(ISTAT,JSTAT, 1 . 
CALL QIO("12,9) 

CALL QIO("12,7) 

PAUSE ' - ADVANCE TAPE TO 

CALL AINIT(ISTAT,IEFM1 ,0) 

TYPE 500, BLINK 

TYPE 7100 

TYPE 500, RESET 

GOTO 901 

END IF 

IF(MIN.EQ. 60) THEM 
CALL ACSTAT(ISTAT,JSTAT, 1 
CALL QIO("12,9) 

CALL QIO{"12,7) 

PAUSE ’ - ADVANCE TAPE TO 
CALL AINIT(ISTAT,IEFM1 ,0) 
TYPE 500, BLINK 
TYPE 7100 
TYPE 500, RESET 


IEFN1 ,0) 

START OF LOW CALS NOW; THEN 
IEFN1 ,0) 

START OF DATA NOW; THEN ' 


GOTO 901 
END IF 

IF (MIN .GE. MAXRUN) GO TO 9100 


i 
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IF PAUSE OR ABORT, END DATA ACQUISITION 


CALL READEF( ij,IEF) 

IF(IEF .GT.O) THEN 
IRSP =IBUF 

IF (IRSP .EQ. I PAWS) THEN 

CALL ACSTAT( ISTAT, JSTAT, 1 , IEFN1 ,0) 

CALL QIO("12,9) 

C CALL QIO("12,7) 

CALL AINIT ( ISTAT, IEFN1 ,0) 

GOTO 7000 
ELSE 

IF (IRSP .EQ. lABORT) GO TO 9100 
END IF 

C SOMETHING ELSE TYPED IN; QUE I/O REQUEST AGAIN 

CALL QI0(IREAD,7, 4 , , lOSB , IPAR , IDSW) 

END IF 
IK=K 

IF(i;<.EQ.1) K=2 
IF(IK.EQ.2) K=1 
IMINS=IMIMS+1 
GOTO 80001 

C CLEAN UP AND EXIT 

9100 CONTINUE 

TYPE 500, RESET 
CLOSE ( UN I T= 2) 

CALL QI0("12,7) 

TYPE 91000 

91000 FORMAT (//////, 10X, 'SELECT ONE OF THE FOLLOWING [ENTER 1 OR 2]' 

1 /20X,'1 TO SAVE DIGITIZED FILE 

2 /20X,'2 TO DELETE DIGITIZED FILE ',$) 

ACCEPT 617,ANS 

IF(ANS.EQ.'I’) THEN 
CLOSE (UNIT=1) 

ELSE IF(ANS.EQ. '2' ) THEN 
CL0SE(UNIT=1 , STATUS= ' DELETE ' ) 

ELSE 

GO TO 9100 
END IF 

CALL ACSTAT( 1ST AT3, JSTAT, 1 , 5 ,0) 

CALL QIO("12,7) 

CALL QIO("12,9) 

END 

SUBROUTINE DIRECT (TCOUNT,K) 

C STORAGE SEQUENCE AS FOLLOWS: 

C FIRST - RAW RESPIRATION 

C SECOND - DELTA Z 

C THIRD - DZ/DT 

C FOURTH - ECG 


34 


c t:<3 for task build 

c 

c Ti(3>FILENAME = FILEMAME 

C TK3>/ 

c 

PARAMETER (IFREQ=200) 
integer buf ( 8*IFREQ) , tcount 
COMMON BUF 

C do 100 tcount = 1,2 

WRITE! 1 ,REC=TCOUMT,ERR-999)( 3UF((K-1 )^J^*IFREQ+ICOUNT) , 

1IC0UNT=1 ,4*IFREQ) 

C TYPE *,'0K OM READ' 

C do 200 i = 1, 200 

C wr i te ( 5 , 1 0 ) i , ( BUF( ICOUNT , 1 ) , IC0UNT= ( i- 1 ) *4+ 1 , i^4 ) 

C200 continue 

C TYPE ^,’OK OM WRITE’ 

100 continue 

10 FORMAT(x,i3,4(2x,i4)) 

20 FORMAT! IX, 15) 

999 CONTINUE 

RETURN 
END 

SUBROUTINE SETFRQ!CFREQ,TFREQ) 

C 

Q»*»»»»3t»»**»*»»*^«*jHt*4f*****^** ************** ■■'********************************** 

C 

C SUBROUTINE GETFRQ!CFREQ,TFREQ) 

C 

C LANGUAGE: FORTRAN-77 

C 

C 

C FUNCTION: 

C This SUBROUTINE will prompt the user for the desired clock frequency to 

C be used in acquiring analog data. 

C 

C 

C OUTPUTS: 

C CFREQ = REAL*4 variable containing the user's desired frequency in hz. 

C TFREQ = REAL*4 variable containing the user's actual frequency in hz. 

C 

C 

C SUBROUTINES REFERENCED: 

C SUBROUTINE CLKFRQ. 

C 

c 

(2 ***»,«******»*»******************************* ********* 

C 

Q »»*»*»jHtX**»»*»»*************»*»**************»***»************ 

c 

C 
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c 

1 


c 


: i ri ^ '-r 


'Tqrrn 


:freo 


IDeclare desired, actual frequencies. 


CONTINUE 

CALL CLKFRQ( CFREQ ,TFREQ) 

IF(TFREO. ME. -999.0) GO TO 2 ’Skip ahead if CFREQ is ok. 


TYPE 9015 

9015 FORMAT( / , IX, ' Bad frequency, please try again’,/) 

GO TO 1 ! Prompt for frequency again. 

C 

2 CONTINUE 

999 RETURN ! Return to caller. 

END 

PROGRAM IMPSTR 

C 

C AN IMPEDANCE CARDIOGRAPH PROGRAM USING DIGITIZED DATA FILES 

C 

PARAMETER (IFREQ=200) 

DIMENSION SL0PE(2) ,ENTRCP(2) 

C COMMON I3?ACE(4*IFRE0 ) , SPACE( 30^IFREQ ) 

REAL*8 DRAW ,IMPACQ , IMPCAL 


C = = = = = = = = = = : = = - = = = = = = = = = - = = = = = = = = = = = = = = = 

C DATA ACQUISITION OF BIOFEEDBACX PARAMETERS 

C ========================================== 


REAL TEMPRT, 

@ ZERO, CFREQ, TFREQ, EPOCH 

INTEGER PREDI A , PREHR , PRETEM , PRESYS , PSTDI A , PSTSYS , NUMDI A , 

@ PSTHR , PSTTEM , RUNMUM , I D A Y , MSN , AGE , SUSC , 

@ TESTIM , MIN , PSTTST , PART , CHANEL , TOTMIN , PRETST , 

@ IMINS,IANS, HIGH( 2,2) ,DIAG(15), 

@ CLEANS(4) ,ICLEAM(2) 

LOGICAL* 1 SEX , ANS , TEMP , DMY , I RSP , I ABORT , I PAWS , TYME AM , STDDEV , 

@ GR0UP(8) ,DATE(7) ,DIRECT(2),0SN(2) ,STRNG(8), 

§ FRMT(12) 

CHARACTER* 1 CMT ( 80 ) , RUN ( 4 ) , RUNT YP , DATF I L ( 1 4 ) , EXP 
CHARACTER*14 CATFIL 
EQUIVALENCE (DATFIL, CATFIL) 

DATA CLEANS/27,91,50,74/ 

C OPEN (UNIT=1 ,NAME='NUDRAW.TEL' ,F0RM='F0RMATTED' ,TYPE='0LD' ) 

C READ(1,2) SLOPE, EMTRCP,C,BSA, CATFIL, MAXREC 

C2 FORMAT (6(G15 . 7) , A 14 , 1 10) 

C CL0SE(UMIT=1) 

9115 CONTINUE 

C TYPE 5020 

5020 FORMAT(/’ EXPERIMENT DESIGNATION (A,B,...): ',$) 

C ACCEPT 5021, EXP 

5021 FORMAT (A1) 

TYPE 5030 

5030 FORMAT(/' GROUP TYPE (maximum of 8 characters): ',$) 

ACCEPT 5031 ,(GROUP(I) ,1=1 ,8) 
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5031 FQRMAT{8A1) 

C TYPE 5040 

5040 F0RMAT(/’ RUM TYPE (maximum of 4 characters): ',$) 

C ACCEPT 504 1 , ( RUM( I ) , I = 1 , 4 ) 

5041 FORMAT (4A1) 

TYPE 5050 

5050 FORMAT(/' RUM NUMBER: ',$) 

ACCEPT *,RUMMUM 
TYPE 5060 

5060 FORMAT(/’ DIRECTION (cc,cw,nd): ’,$) 

ACCEPT 506l,(DIRECT(I),I=l,2) 

506 1 FORMAT (2A1) 

C TYPE 5070 

5070 FORMAT(/' DAY NUMBER: ',$) 

C ACCEPT *,IDAY 

C TYPE 5080 

5080 FORMAT(/' MEW SUBJECT NUMBER: ',$) 

C ACCEPT *,NSM 

TYPE 5090 

500C FORMAT(/' SUBJECT' 'S INITIALS (first and last): ’,$) 

ACCEPT 5091 ,(OSM(I) ,1=1 ,2) 

5091 FORMAT (2A1) 

TYPE 5100 

5100 FORMAT(/’ AGE: ',$) 

ACCEPT »,AGE 
TYPE 5110 

5110 FORMAT(/' SUSCEPTIBILITY (1,2,3): ',$) 

ACCEPT »,SUSC 
TYPE 5120 

5120 FORMAT(/' SEX (m/f ) : ',$) 

ACCEPT 5121, SEX 

5121 FORMAT (A1) 

TYPE 5130 

5130 FORMAT(/' TEST DATE (mmmddyy): ',$) 

ACCEPT 5131 , (DATE( I ) , 1=1 ,7) 

5131 FORMAT (7A1) 

TYPE 5140 

5140 FORMAT(/' TEST TIME (military in hrs and mins, e.g. for 1:30 p.m. 
@ ,' enter 1330): ',$) 

ACCEPT »,TESTIM 
TYPE 5150 

5150 FORMAT(/’ PRE HEART RATE (in beats/min): ',$) 

ACCEPT *,PREHR 
TYPE 5160 

5160 FORMAT(/’ PRE TEMPERATURE: ',$) 

ACCEPT *,TEMPRT 

PRETEM = IIFIX(TEMPRT*10.0) 

TYPE 5170 

5170 FORMAT(/' PREB.P. (sys,dias ,e.g. 120,80): ',$) 

ACCEPT *,PRESYS,PREDIAS 
TYPE 5180 

5180 FORMAT(/' PRE/POST TEST BASELINE (in minutes): ',$) 

ACCEPT *,PRETST 
PSTTST = PRETST 
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CATFIU r : 14) = ’ . IMP’ 

CALL STORE ( GATE IL , IFLAG , GROUP , RUNMUM, DIRECT ,OSN, ACE ,SUSC , SEX , 

1 DATE , TESTIM , PREHR , TEMPRT , PRETEM, PRESYS , PREDI A , PRETST , PSTTST ) 
TYPE 500, CLEANS 
500 F0RMAT(X,4A1) 

END 

SUBROUTINE STORE (CATFIL, I FLAG, GROUP, RUNNUM, DIRECT, OSN, AGE, 

1 SUSC , SEX , DATE , TEST I M , PREHR , TEMPRT , PRETEM , PRESYS , PREDI A , PRETST , 

2 PSTTST) 

C 

C AM IMPEDANCE CARDIOGRAPH PROGRAM USING DIGITIZED DATA FILES 

C 

DIMENSION SL0PE(2) ,ENTRCP(2) 

COMMON DATCAR(4,2) 

REAL*8 DRAW 

C = = = = = = = = r = = = = r = = = = = = = = = = = = = = = = = = = = = = 

C DATA ACQUISITION OF BIOFEEDBACK PARAMETERS 

C : = = ^ = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =: = = = 

REAL TEMPRT, MEAM(4) ,SD(4) , 

@ ZERO, CFREQ,TFREQ, EPOCH 

INTEGER PRED I A , P REHR , PRETEM , PRESYS , PSTD I A , PSTS Y3 , NUMD I A , 

@ PSTHR , PSTTEM , RUNNUM , I D AY , NSN , AGE , SUSC , 

§ TESTIM , MIN , PSTTST , P ART , CHANEL , TOTMIN , PRETST , 

§ IMINS,IANS, HIGH( 2,2) ,DIAG(15), 

§ CLEANS(4),ICLEAN(2) 

LOG I C AL » 1 SE X , ANS , TEMP , DMY , I RSP , I ABORT , I P AWS , TYME AN , STDDE V , 

§ GR0UP(8) ,DATE(7) ,DIRECT(2) ,0SN(2) ,STRNG(8) , 

@ FRMT(12) 

CHARACTER* 1 CMT ( 80 ) , RUNTYP , DATFIL( 1 4 ) , EXP 
CHARACTER*4 RUN 
CHARACTER*14 CATFIL , IMPFIL 
EQUIVALENCE (DATCAR,MEAN) , (DATCAR( 1 ,2) ,SD) 

DATA DRAW/’DRAWRMV , ICLEAN/27 , 99/ 

DATA cleans/27,91,50,74/ 

C IFLAG=0 

C TYPE 500, CLEANS 

500 F0RMAT(X,4A1 ) 

OPEN ( UNIT= 1 , STATUSr ’ OLD ' , NAME= ' DW 1 : NUDRAW . TEL ' , FORM= ' FORMATTED ' ) 
READ( 1,2) SLOPE , ENTRCP , C , BSA , IMPF I L , MAXREC 
2 F0RMAT(6( G15 . 7) , A14 , 1 10) 

CL0SE(UNIT=1 ,DISPOSE='SAVE' ) 

IMPFIL(12:14)='IMP' 

MINrMAXREC/4 ! THERE ARE 4 EPOCHS PER MINUTE 
N=1 .'SET FOR FIRST DISKETTE 
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c 


ASSEMBLE DATAFILE fJAME 


continue 

CALL ERRSET(63, .TRUE. , .FALSE. , .FALSE. , .FALSE. ,MA.X) 
CALL ERRSET(36, .TRUE. , .FALSE. , .TRUE. , .FALSE. ,MAX) 
CALL ERRSET ( 29 TRUE FALSE . TRUE . , . FALSE . , MAX ) 
EX?=IMPFIL(5:5) 

DEC0DE(2,9112,IMPFIL(6:7)) WSN 
DECODE(2,9112,IMPFIL(8:9)) IDAY 


IF 

{ IMPFIL( 10:10) 

.EQ. 

'O' ) 

RUN( 1 :4)='BSLN' 

IF 

(IMPFIL{ 10:10) 

.EQ. 

’ 1') 

RUN( 1 :4) = 'CSSr 

IF 

(IMPFIL( 10: 10) 

.EQ. 

'2' ) 

RUN( 1 :4)='TRAN' 

IF 

(IMPFIL( 10:10) 

.EQ. 

'3') 

RUN(1:4)=' ' 

IF 

(IMPFIL( 10:10) 

.EQ. 

'4' ) 

RUN( 1 :4)='TASK' 

IF 

(IMPFILC 10:10) 

.EQ. 

'5') 

RUN( 1 :4)='AMBL' 

IF 

(IMPFILC 10: 10) 

.EQ. 

'6') 

RUM(1:4)=' 

IF 

(IMPFIL( 10:10) 

.EQ. 

'7' ) 

RUN( 1 :4) = 'VARD' 

IF 

(IMPFIL( 10:10) 

.EQ. 

'8' ) 

RUN( 1 :4)='DRUM' 

IF 

(IMPFIL( 10:10) 

.EQ. 

'9') 

RUN{ 1 :4) = 'ZERO' 


C.ATFIL(1:1) = 'D' 

CATFIL(2;2) = 'Z' 

CATFIL(3:3) = T 
CATFIL(4:4) 

CATFIL(5:10)=IMPFIL(5:10) 

9111 FORMAT(A2) 

9112 FORMAT( 12) 

CATFIL( 11:11) = ' 

C 

C OPEN DATA FILES 

C 

TYPE 500, CLEANS 
40 TYPE 9113 

9113 FORMAT(//////////////////' PLEASE INSERT DATA FLOPPY WITH AT 

1 LEAST 60 CONTIGUOUS FREE BLOCKS'/' INTO DRIVES DZ1 & DZ2.', 

1 ' WHEN READY, PRESS RETURN',/////////) 

ACCEPT 10,ITEMP 
10 F0RMAT(I5) 

C 

9510 CONTINUE 

C 

CATFIL( 12: 12) = '!' 

CATFIL{13:13) = 'M' 

CATFIL( 14:14) = 'P' 

IRECrO 
CALL ERRSNS 

OPEN (UNIT=1, TYPEr'OLD', NAME=CATFIL, FORM =' FORMATTED ' , 

1 ERR=45 ) 

GOTO 100 
45 CONTINUE 
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O O Ci 


I, linaL .'! 


OPEN 

rwT-Uf 


(UNIT: 


a-: -:. - : ■ .y- ' i 

IZP.'/EM) 

o 'j “"p E-J nr \j 

''tY?e'=^MEW’ , NAMErCATFIL, FORM= ' FORMATTED’ , 

) 


GoTo 40 

END IF 

ICC CONTINUE 


STORE KEY FIELD INFO TO DATA FILE 


TYPE *, 'EPOCHS, RUN DURATION (MINUTES)= ',IMINS,MIN 
WRITE( i ,9120) EXP,(GROUP(I) ,1=1 ,8),(RUN ), 

1 RUNNUM,(DIRECT(I),I=1 , 2 ) , IDAY ,NSN, (OSN( I ) , 1= 1 ,2) , 

2 AGE , SUSC , SEX , ( DATE( I ) , I = 1 , 7 ) , TESTIM , PREHR, PRETEM, PRESYS , PREDI A , 

7 oQ/T^uo ^ PSTTEM , PST3YS , PSTD I A , MIM , PRETST , 

4 MIN-PRETST-PSTTST, PSTTST.NUMDIA, 1 , 8 ,DIAG(1) 


9 i 20 


FORMATi Ai .X,8A1 ,X, A4 , X , 12 , X , 2A1 , 2( X , 12 ) , X , 2A1 ,X , 12 , X , 1 1 , 

1 X, AI ,X,7A1 ,X. I4,X, I3,X, I4,X, I3,X, 13, 

2 X , 1 3 , X , 1 4 , 5 ( X , 1 3 ) . 5 { X , 1 2 ) ) 


DO 9130 I = 2,14 

WRITE( 1 ,9140) DIAG(I) 

9130 CONTINUE 

9140 F0RMAT(109X,I2) 

WRITE( 1 ,9150) DIAG( 15) ,MAXREC 
9150 FORMAT(109X,I2, 13) 


FRMT(2) = 'O’ 

FRMT(3) = '6' 

ZERO =0.0 
IERR=0 

9650 CONTINUE 

IREC=IREC+1 

IF(IREC.GT.MAXREC) GOTO 9165 
D TYPE IREC= ',IREC,' MAXREC= ' ,MAXREC 

CALL DRECT ( IREC, lERR, IMPFIL) 

IF(IERR.EQ.I) GOTO 9160 

WRITE( 1 ,9550) MEAN,0. ,0.,0.,0.,SD,0.,0.,0.,0. 

C WRITE( 1 ,9550) (DATCAR(I , 1 ) , 1=1 ,4) ,0. ,0. ,0. ,0. , 

C 1 {DATCAR(I,2),I=1,4),0.,0.,0.,0. 

9550 F0RMAT(8(X,F9.4)) 

GOTO 9650 
9160 CONTINUE 

WRITE! 1 ,9550) 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 ., 0 . 
GOTO 9650 
9165 CONTINUE 

CL0SE(UNIT=1) 

CATFIL(12:12) = 'C 
C DATFIL( 13)=’ A' 
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om\n \t I ’ 

OF POOR QUALITY 


C 

c 

C 9600 

C9700 

C 


1 

1 


DA'^riLC 1 4 ) = ' L ’ 


OPEM (UMIT=;, 
ERR=40 


TY?E='MEW’ , NAME-CATF IL, FORMz ' FORMATTED ' 


DO 9700 r=;,2 

WRITE( 1 ,9600) HIGH( 1 . 2 ) , HIGH( 1 , 1 ) ,CALVAL( 1 , 1 ) ,CALVAL(2, I ) 
FORMAK 2(X, 1 10) ,2(X,F10.4)) 


CONTINUE 
CL0SE(UNIT=1 ) 
IF(IFLAG.EQ.O) GOTO 9875 


CATFIL( 13:13)='0' 

CATFIL( 14:14)=’M' 

OPEN (UNTT=1, TYPE: 'MEW, 
ERR=40 ) 


NAME 


rr 


, FORM: ' FORMATTED' , 


OPEN (UNIT:2, TYPEWOLD', MAME= 'SY: TEMP. COM' , FORM: ' FORMATTED 
READONLY ) 

REWIND 2 


J:0 

KCMT:30 
9750 CONTINUE 
J:J-'KCMT 

IF( J .GT . KCHAR ) KCMT:KCHAR- J+KCNT 
READ(2,9800) (CMT( I ) , 1:1 ,KCMT) 

WRITE! 1 ,9801 )(CMT( I ) , 1=1 ,KCMT) 

9 SOO FORMAT! 30A I ) 

9801 FORMAT! X,80A1) 

IF!J.GE. KCHAR) GOTO 9850 
GOTO 9750 
9850 CONTINUE 

CL0SE(UNIT=1) 

CLOSE! UNIT=2) 

9875 CONTINUE 

IF!N.EQ.2) GOTO 99990 
N:2 

CATFIL!3:3)='2' 

GOTO 9510 
99990 CONTINUE 

TYPE ARE YOU FINISHED WITH THE DIGITIZED DATA ?[Y/N] ' 
ACCEPT 10000, ANS 
10000 FORMAT! A1) 

IF!ANS.EQ. 'Y' ) THEM 

TYPE THEM YOU WANT TO DELETE THE DIGITIZED FILE ?[Y/N]’ 
ACCEPT 10000, AMS 
IF!ANS.EQ. 'Y' ) THEN 

OPEN !UMIT:2,NAME:IMPFIL ,rYPE:'OLD') 

CLOSE !UMIT:2,DISP0SE: 'DELETE' ) 

IMPFIL! 11:14):' .DIG' 

OPEN !UNIT:2,NAME:IMPFIL , TYPE: 'OLD') 

CLOSE !UMIT:2, DISPOSE: 'DELETE' ) 

END IF 
END IF 

IF!IFLAG.EQ. 1) THEN 

OPEN !UNIT:2,NAME:'SY:TEMP.C0M' ,TYPE='OLD' ) 

CLOSE !UNIT=2,DISP0SE:'DELETE' ) 

END IF 
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c 

c 

gnoo 

99999 


c 

c 

c 

c 

c 


OFEM (UM:T:1,G7ArU3='CLD’ ,?IAME=’D'vi 
CL03E(lJW: r^ i ,D[S?OSE=' DELETE' ) 


TYPE 91100 

FORMAT (/' SUCCESSFUL COMPLETICM ') 

COMTIMUE 

RETURN 

END 

SUBROUTINE DRECT (TCOUNT, lERR, IMPFIL) 
STORAGE SEQUENCE AS FOLLOWS: 

FIRST - RAW RESPIRATION 
SECOND - DELTA Z 
THIRD - DZ/DT 
FOURTH - ECG 


;RM- ' FOkMA I in.:. 


c TKB FOR TASK BUILD 

C 

c TKB>FILENAME-FILENAME 

C TKB>MAXBUF=1600 

C TK3>// 

c 


PARAMETER (IFREQ^EOO) 

C integer buf (4^IFREQ, 1 ) 

INTEGER tcount 

C DIMENSION X(15*IFREQ),Y{15*IFREQ) 

DIMENSION DATCAR(4,2) 

COMMON DATCAR 

C EQUIVALENCE ( DATCAR, BUF) 

CHARACTER*14 IMPFIL 


OPEN (UMIT=2, status='old' , MAME=IMPFIL , 

1 FORM: ' UNFORMATTED ' , access: ' direct recl=8 ) 


C 

c 

c 

c 

C200 

C 

100 

990 

999 

10 

20 


do 100 tcount : 1,2 

read(2,REC:TCOUNT,ERR:990)((DATCAR(I,J),I=1 ,4),J=1 ,2) 

TYPE *,'0K ON READ' 
do 200 i = 1, 200 

write(5, 10)i,(BUF(IC0UNT, 1),IC0UNT=(i-1)»4+1 ,i*4) 

continue 

TYPE *, 'OK ON WRITE' 

continue 
GOTO 999 
CONTINUE 
IERR=1 

CL0SE(UNIT:2, STATUS: ' SAVE ' ) 

FORMAT(x,i3,4(2x,i4)) 

FORMAT( IX, 15) 


RETURN 

END 

PROGRAM DRAWED 

PARAMETER (IFREQ:200) ! SET SAMPLING SIZE FOR DIGITIZED DATA 
DIMENSION ISTAT(2) ,IRES(75) ,OUTPUT( 10,70), 
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o n o o o 


1 


V y 

c 

4000 

4500 

1 


SLC?E(2) ,i:'JTRCP(2) 

CHARACTER*! AMS 

CHARACTER* 1 ICHRPK( 2 ) , I BUF , ISTOP , IRE PET 
CHARACTER*4 EPOCH ,ETIME 
VIRTUAL IDATA(9000) 

DIMENSION Y( 1450) ,X( 1000) ,NCHAN(3) ,SMIN{3) ,SMAX(3) ,DATCAR(5,50) , 
DATSD(4),IPAR(6) 

INTEGER*2 IRESP(IFREQ) 

CHARACTER* 14 CATFIL 

INTEGER LPTS( 3 ) , I PEAKS( 50 ) , IOSB( 2 ) , CLEANS! 2 ) 

BYTE 1233. 173, 1162, IH0ME{4) 


CHARACTER*2 IPT,IPB 


COMMON /PLOT/X,Y 

EQUIVALENCE (DATCAR,Y( 1201 ) ) ,(Y(501 ) .OUTPUT) 

DATA NCHAM/3,4, 1/,SMIN/350. , 150. ,0. / ,SMAX/550. , 350. , 150 . / 

DATA LPTS/3 ,3 ,10 / , MSMFRQ/ 100/ ,NREC/ 15/ , IREAD/"001000/ 


DATA A/.01/,B/0./,C/0./,IBRAC/' ]'/,ICHAN/3/ 
DATA IGATE/0/ , IM/0/ , I233/"233/ , I73/"73/ , 1 162/ 


" 162 / 


DATA IHOME/155, 63, 54, 108/, CLEANS/27. 99/ 

CALL GETADR(I?AR(1),I3UF) !GET START ADDR OF I3UF & STOP. 
IPAR{2)=1 !SET TO BYTE SIZE OF IBUF 
TYPE 4000, CLEANS 
F0RMAT(X,4A1 ) 


IPAR( 1 ) 


TYPE 4500 

FORMAT ( 10X,60( '*'),///// ,22X, ' DATA REDUCTION 4 EDITING , 

/////, 10X,60('*’),////) 


READ FILE CONTAINING SLOPES, INTERCEPTS, C AND BSA 
WHERE C = RH0*L**2 

BSA = BODY SURFACE AREA 


OPEN {UNIT=1 ,NAME='DW1 :NUDRAW.TEL' , FORM= ' FORMATTED ' ,TYPE=’OLD' ) 
READ! 1,2) SLOPE, ENTRCP,C, BSA, CATFIL, MAXREC 
2 F0RMAT!6! G15.7),A14,I10) 

CLOSE !UMIT=1) 

IDZDT0=-ENTRCP!2)/SL0PE!2) ! GET AD VALUE FOR DZDT=0 
IDZDT1=!1.-ENTRCP!2))/SL0PE!2)! GET AD VALUE FOR DZDT=1 
100 CONTINUE 

TYPE *,’ STARTING TIME CODE FOR DATA ON TAPE !eg. enter 300 for', 
1 ' 3 : 00 )’ 

ACCEPT *,ISTREC 

TYPE *, ’ EPOCH FOR START OF DATA REDUCTION ' 

ACCEPT *,ISTEPC 

TYPE *, ’ ENTER LAST DATA EPOCH ' 

ACCEPT *,IENEPC 

c 

C KREC SETS THE STARTING RECORD FOR DATA IN DIGITIZED FILE 

C INDREC SETS THE ENDING RECORD FOR DATA IN DIGITIZED FILE 

C IDKREC SETS THE 15 SECOND EPOCH NUMBER ! DISPLAYED TO SCREEN) 

C AND USED TO STORE RESULTS IN 15 SEC EPOCHS 

C 
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L'i£r.£C=(I£M£?C)»{5>6l 
MAaHc.C = { INDREC-6 1 ) / 1 5 
I OK P£C = I STS ?•' 

CArFIL{ 11:14)=’ . DIG ' 


105 
C 
1 0 


c 

c 


c 

125 

0998 

C 

C 

C 

c 

c 

c 




CONTIWUE 

CALL QI0(IREAD,7,4, , I0S3 , 1 PAR , IDSW) 

CONTINUE 

XMIM=0. 


A.'! HA 


= iirREQ “NREC 


:ALL CGL( 90) UNITIALIZE GRAPHICS 


CALL 
CALL 
CALL 
CALL 
CALL 

'"'.r'crv /r^ATT' 


CGL( 92) !NEW FRAME 

CGL( 103, 'DW1 :FILE2.GID' , 13) ! INITIALIZE PLOT FILE 
CGL( 105, 'DW1 :FILE2.GID' , 13) fSELECT PLOT FILE 

CGL( 80,XMIM,XMAX,0. ,600. ) !SET WINDOW TO DEFAULT VALUES 

CGL(86,0) .'SET ORIGIN 




^,T _ r* 0 'V rr T T 


FOFMr ' UNFORMATTED ' 

P 192: ^ ^ ' T 17 IT r 

IR£C = KR£C-hI- 1 


!0?EM DIGITIZED DATA 


rr p 


access^ ’direct ' , reci=i|^IFREQ) 


read{ P REC= IR aC , ERR=999 ) ( ( IDATA( J+( 1-1 )^IFREQ) , ID, ID, ID) , 
1 U- 1 , xtREQ) 

CONTINUE 

CL0SE(UNIT=1 , STATUS= 'SAVE' ) 

M=0 

NPTS=LPTS(3) 

DO 135 L=1,IFREQ^NREC,NPTS 

M=M+1 

X(M)=L 

Y(M)=IDATA(L) 

CuNT J.NUE 

M=0 !SET COUNTER OF RESPIRATION VALUES 
IFRQRC=IFREQ^NREC !SET AREA SIZE FOR DIG ( IDATA)BLOCKS 
DO 150 I=1,MREC 
IREC = r.:REC-^I-1 


STORAGE SEQUENCE OF DIGITIZED DATA AS FOLLOWS: 

FIRST - RAW RESPIRATION - STORED IN RESP FROM 1 TO IFREQ 
SECOND - DELTA Z - STORED IN IDATA FROM 1 TO IFRQRC 
THIRD - DZ/DT - STORED IN IDATA FROM UIFRQRC TO 2*IFRQRC 
FOURTH - ECG - STORED IN IDATA FROM 1+2*IFRQRC TO 3*IFRQRC 

read( 1 ,REC=IREC ,ERR=999) ( (IRESP( J) , 

1 IDATA( J+(I-1 )»IFREQ) , IDATA( IFRQRC+J+( 1-1 )»IFREQ) , 

2 IDATA(2»IFRQRC+ J+( I- 1 ) *IFREQ) ) , J= 1 , IFREQ) 

PRINT *,(IECG(J+(I-1)*IFREQ ) , J= 1 , IFREQ) , J , I 

STORE A SUBSET( EVERY LPTS(3) VALUE) OF RESP IN Y FOR PLOTTING 


DO 175 J=1, IFREQ, LPTS(3) 
M=M+1 


ORIGINAL 

OF POOR Oti.AL^rY 
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150 




COMTIN'JZ 

COMTIMUE 

CLOSZ' uMIT- V. STATUS: 'SAVE' ) [CLOSE DIG. DATA FILE 


C 


PLOT FIRST n£SPIRATI0N(K=3) THEN EC(:(K=2) AND LAST DZ/DT(K=1 


DO 3000 K=3J,-1 

Y3MIM=3MIN(r:) [PARTITION SCREEN ARE., MIN FOR K. PLOT 
YSMAX=SMAX(K) [PARTITION SCREEN AREn MAX FOR K PLOT 
MPT3=LPTS(K) 


C SET PLOT ARRAY PAIRS (X,Y) 

C 

IF(K.EQ.3) THEN 

DO 275 L=1,M [SETUP X VALUES FOR RESPIRATION PLOTTING 
X(L):1+(L-1 )^NPTS 

CL 


DO 300 L:1.IFREQ ^NREC,N?T3 [SETUP (X,Y) PAIRS FOR ECG AND DZ/DT PLO 

M:M+1 

X(M)rL 

Y ( M ) : I DAT A ( I FRQRC ^K+L ) 

300 CONTINUE 
END IF 

CALL INDEX(XMIN,XMAX,YMIN,YMAX,M ) [DETERMINE PLOT LIMITS 
IF(K.EQ.I) THEN 

IF( IDZDT1-IDZDT0.lt. (YMAX-YMIN)/3) THEN [REDONE LIMITS IF NEEDED 
YMAX = I DZDT 1 + 1 DZDT 1 - 1 DZDTO 
YMIN=IDZDT0-(IDZDT1-IDZDT0) 

END IF 

IF(IDZDT1 .GT.YMAX) YMAX = I DZDT 1 

DZDTO= ( IDZDTO-YMIN) /( YMAX-YMIN) *( YSMAX-YSMIN)+YSMIN [ SCALE DZDTO 
DZDTU(IDZDTT-YMIN)/(YMAX-YMIN)*(YSMAX-YSMIN)+YSMIN [SCALE DZDT1 

PLOT DZ/DT CALIBRATION VALUES FOR OESERVATION 

CALL CGL( 1,X(1), DZDTO) [MOVE "PEN" TO POSITION 

CALL CGL( 4,XMAX, DZDTO) [DRAW TO POSITION 

CALL CGL( 1,X(1), DZDT1) [MOVE "PEN" TO POSITION 

CALL CGL{ 4,XMAX, DZDT1) [DRAW TO POSITION 

END IF 

DO 322 L=1,10 

TYPE *,X(L),Y(L) ,YMIN,YMAX,L 

Y(L)=(Y(L)-YMIN)/(YMAX-YMIN)*( YSMAX-YSMIN)+YSMIN 
TYPE *,X(L) ,Y(L) ,YMIN,YMAX,L 
C322 CONTINUE 

DO 325 L=1,M [SCALE Y 
C PRI.NT *,X(L) ,Y(L) ,YMIN,YMAX,L 

Y ( L ) r ( Y ( L ) - YM I N ) / ( YM AX - YM I N ) ^ ( YSMAX - YSM I N ) + YSMIN 
C TYPE *,X(L) ,Y(L) ,YMIN,YMAX,L 

325 CONTINUE 

C FACT=5. 732767. 


Of. 
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C XOLD = XUOO)+ ■ ./2^.0. 

c m'/::MREC":freq 

MYPLT: IFRCEC /MFTS 

CALL CGL^F.X. Y.MYFLT) ! PLOT MYPLT X.Y PAIRS 
IF (X.EQ.i) THEN ! PLOT EPOCH NUMBER AND START TIME OF EPOCH 
IEPOCH=IMOD( KREC-6 1 , 60 )+( KREC-6 1 ) /60^ 100+ISTREC 
ENCODE (4,32500,ETIME)IEP0CH 
ENCODE ( 4 , 32500 , EPOCH) IDKREC 
32500 F0RMAT(I4) 

CALL CGL ( 1 ,XMIM.SMIM( 1 ) } 

CALL CGL ( 1 5. EPOCH, 4) 

CALL CCL { 1 ,XM!M^240. ,SMIM( 1 ) ) 

CALL CGL (16,ETIME,L) 

END IF 

500 CONTINUE 

3000 CONTINUE !END OF PLOTTING 
C 

C NEXT CALL PEAK DETECTION 

n 

D OPEN (iiNIT = 2,STATUS='MEW , NAME= ' DW1 : IMPTEMP .TST ' , 

D 1 FORM= ’ FORMATTED ACCESS SEQUENTIAL ' ) 

CALL VOL( IFRORC .NPEAK3, IDA.TA, lADMIN, IGATE, IM, IDZDTO, IDZDn ) 

C IDZDT - INPUT OF DZ/DT DATA 

C IFRQRC - NUMBER OF DZ/DT DATA POINTS 

C OUTPUT! 3, I) - ARRAY CONTAINING TIME VALUES FOR DETECTED PEAKS 

C NPEAKS - NUMBER OF PEAKS DETECTED 

IF(NPEAKS.GT.O)THEN 

CALL CGL( 12,4,0,0) !SET LINE STYLE TO DOTTED 
C 

C PLOT PEAK INDICATORS 

C 

DO 630 1=1, NPEAKS 
YLOW=SMIN(3)+40 
XL=(0UTPUT(3,I)-1 . ) 

CALL CGL( 1 ,XL,SMAX( 1 ) ) ! POSITION TO X=XL ,Y=SMAX( 1 ) 

CALL CGL(4, XL,YLOW) IDRAW LINE TO X=XL,Y=YLOW 
ENCODE (2,610,ICHRPK),I 
610 FORMAT! 12) 

XL0W=XL-10. 

DO 620 J=1,2 

CALL CGL! 1 ,XLOW,YLOW ) 

CALL CGL! 1o,ICHRPK!J) , 1 ) IWRITE PEAK NUMBER J 
YL0W=YL0W-20. 

620 CONTINUE 

680 CONTINUE 

CALL CGL! 12, 1 ,0,0) IRESET LINE STLYE TO SOLID 
C 

C USE A SUBROUTINE HERE TO CALCULATE THE TIME WHEN 

C dZ/clT=0 BEFORE EACH PEAK ALSO CALCULATE STROKE VOLUME 

C 

CALL TFACT!NPEAKS,IFREQ ,C,NY ,K,TYMEAN, 

1 SLOPE, ENTRCP,BSA,IDATA,KREC, IDZDTO, I ADMIN, CATFIL.IDKREC) 

C ELSE 


Of 
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c 

c 

c 

c 

c 

c 

c 

C700 

D 

C 

C 

c 


40000 


725 

C 

C 

C 

c 

c 

c 

c 


SET MEAM AMD SD TO ZERO 
END IF 

IF( ICOUMT.GT. 1 ) THEN 
CALL CGL( 12,4,0,0) 

DO 700 I=1,IC0UNT-1 
XL=600. *1/200. 

CALL CGL( 1 ,XL,YMAX ) 

CALL CGL(4,XL ,YMIN ) 
CONTINUE 
CL0SE{UNIT=2) 

SET UP ONE LINE DIALOG AREA NOW 


L20=49 

L24=49 

WRITE(5,40000) I233,L20, I73.L24, 1 162 
FORMAT( 1X,5A1 ) 

DO 725 1=1,50 IZERO ALL IPEAKS - ARRAY OF USER 

IPEAKS( I )=0 

CONTINUE 


SELECTED 


PEAKS 


QUIRY EVENT FLAG 4 TO SEE IF USER REQUEST MADE 


CALL READEF(4,IEF) 

IF(IEF.GT.O) THEM ! REQUEST MADE VIA KEYBOARD 
END IF 

CALL QIO("12,7) ! CANCEL QUE I/O REQUEST 


DIALOG SECTION 


ONE 

PLOT FILE 
'.TERMINATE PLOT FILE 


TYPE SKIP THIS EPOCH ? [ENTER Y/N]' 

ACCEPT 10000, ANS 
IFCANS.EQ. ' Y’ ) THEN 

KREC=KREC+MREC ! INCREMENT RECORD INDEX 
MAXREC=MAXREC-1 ! REDUCE MAX RECORD SIZE BY 
CALL CGL(106,'DW1:FILE2.GID' ,13) 'DESELECT 
CALL CGL( 104, 'DW1 :FILE2.GID' ,13) 

CLOSE (UMIT=1,DISP0SE=' DELETE') 

IF (KREC.GE.INDREC) GOTO 99999 
GOTO 105 
END IF 

TYPE REDO ANALYSIS ? [ENTER Y/N]' 

ACCEPT 10000, ANS 
IF(ANS.EQ. ) THEN 

TYPE ENTER GATE FACTOR, MINIMUM DIFFERENCE FACTOR NOW 

ACCEPT *,IGATE,IM 

CALL CGL(106,'DW1:FILE2.GID',13) 

CALL CGL(104,'DW1:FILE2.GID',13) 

CLOSE (UMIT=1 ,DISPOSE=' DELETE' ) 

GOTO 10 
END IF 

TYPE OK TO CALCULATE HEATHER 

ACCEPT 10000, ANS 


INDEX ? [Y/N]' 
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7<i50 


7025 


CALL HICALCvMPEAKS, I FREQ 
ELSE 

DO 7250 r=!,.'J?EAKS 
DATCAR(U. l)=o.o 


, I DATA 


,KREC,NREC,CATFIL) 


CONTIi'iu'E 
END IF 

CALL CGL( 106, 'DW 1 :FILE2.GID' , 13 ) 

CALL CGL( 104, 'DW1 ;FILE2.GID' , 13 ) 

CALL CGL(93,INAM,IC0DE) ! REPORT ERROR 
It ( xCODE.NE.O) THEM 

TiPE IMAMx IC0DE= ' , ICODE 

CALL CGL( 91 ) ! TERMINATE GRAPHICS 
END IF 
CONTINUE 


C 


7725 

730 

740 

r 

C 


750 

775 

C 


C 


810 


820 


CALL ERRSNS 

PRINT NPEAKS= ' ,NPEAKS 

* ' ( - AKS ( 1 1 ) , 1 1 = 1 , MPEAKS ) 

IiPt LIST SELECTED PEAKS MOW (SEPERATE W/ COMMA )' 
READ(5.7725'' ( IPEAKSfl ), Iz 1 , MPEAKS) 

FOnMAT( <NPEAKS > 13 ) 

CONTINUE 

CONTINUE 


PRINT MFEAKSz ' ,NPEAKS 
PRINT *, (IPEAKS(II),II=1 ,NPEAKS) 
KPEAKS =0 


DO 750 1 = 1, MPEAKS .'DETERMINE NUMBER OF PEAKS SELECTED 

IF(IPEAKSd) .EQ.O) GO TO 775 

KPEAKS=KPEAKS +1 

CONTINUE 

CONTINUE 

TYPE *, - KPEAKS= ',KPEAKS , (IPEAKS( I ) , 1= 1 ,KPEAKS) 
IF(KPEAKS.GT.O) THEN 


PRINT *,IDKREC,( IPEAKS( I ) , ( DATCAR(K, IPEAKS( I) ) ,K=1 ,4) , 1 = 1 KPEAKS) 
CALL STRVOL(KPEAKS,DATSD,IPEAKS) ’ I, SPEAKS) 

VALUE, CARDIAC OUTPUT, CARDIAC INDEX, HEATHER INDEX 

UU o I 0 I = 1 , u 

DATCARd, 1)=0.0 ! MEANS 
DATSD(I)=0.0 ! STANDARD DEVIATIONS 
CONTINUE 
END IF 


CALL ERRSNS(IERR) 
IF(IERR.NE.O) GOTO 7025 
CONTINUE 


WRITE RESULTS OF EPOCH( =IDKREC) CALCS TO INTERNAL FILE 

CATFIL(11:14) = MMP' 

OPEN (UMITrl ,STATUS= 'OLD' , MAME=CATFIL , 

1 F0RM= ' UNFORMATTED ' ,ACCESS=' DIRECT ',RECL =8 ) 

WRITE( 1 , REC=IDKREC ) ( ( DATCAR( K , 1 ) , K= 1 ,4 ) ) , DATSD 
CATFIL( 11:14)=' .DIG' 

CL0SE(UNIT=1) 

CL0SE(UNIT=2) 
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10000 


99999 


C 

C 


C 


C 


300 

C 

C 


C 


c 


TYPE TO YOU WAMT A COPY ? [Y/M]' 

ACCEPT 10000, AW3 
F0RMAT(A1 ) 

IF (AMS.EQ.’Y’) THEM 
IRES( l)=i 

CALL CPRMTlISTAT, IRES, 'DW1 :FILE2.GL'' , 13) 

ELSE 

OPEN (UNIT=1 ,STATUS='OLD' ,NAME='DW1 :FILE2.GID' ) 
CLOSE (UNIT= 1 , DISPOSE= ' DELETE ' ) 

END IF 


TYPE ' EXIT ? [Y/M] ' 

ACCEPT 1 0000 , AMS 
IF(AMS.EQ. ’ Y' ) GOTO 99999 

IDKREC = IDKREC+1 ! INDEX EPOCH COUNTEF- FOR PLOT AND RESULTS FILE 
KREC=KREC+NREC ! INDEX DIGITIZED DATA FILE RECORD COUNTER 
IF(KREC.GE. IiMDREC) GOTO 99999 !AT END OF DATA FILE YET 
GOTO 105 
CONTINUE 


TYPE 4000, CLEANS 

OPEN (UNIT-1 ,NAME='DWi :NUDRAW.TEL' , FOR.M^ ’ FORMATTED ' ,TYPE=’OLD' ) 
WRITE(1,2) SLOPE, EMTRCP,C,BSA,CATFIL,MAXREC 'RECORD MAXREC VALUE 
CLOSE (UMITrl) 


WRITE(5, 40000) 1233,173,1162 

TYPt, EXST CALLED ISTAT= ' , IS TAT 

CALL EXST (I ST AT) .'SET FLAG FOR EXIT TO PARENT TASK 

STOP 

END 


SUBROUTINE INDEX (XMIM,XMAX,YMIN,YMAX,NPTS) 

PARAMETER (IFREQ=200) 

DIMENSION X( 1000), Y( 1450) 

COMMON /PLOT/X,Y 

PRINT *,K,MPTS,(Y(I),I=1,NPTS) 

SX=639./MPTS 
YMIN=Y( 1) 

YMAX=YMIM 
DO 300 1=1 ,MPTS 
PRINT *,Y(I),I 
IF(Y(I).LT.YMIN) YMIN=Y(I) 

IF{Y(I).GT.YMAX) YMAX=Y(I) 

CONTINUE 

SY=479./{IYMAX-IYMIN) 

ENY=-479.*IYMIN/(IYMAX-IYMIN) 

XMIN=0. 

XMAX=X(NPTS ) 

PRINT *,Y,YMIN, YMAX,NPTS 

RETURN 

END 

SUBROUTINE VOL(NY,KPEAKS, lY, lADMIN, IGATE, IM, IDZDT0,IDZDT1 ) 
DO PEAK DETECTION USING DIGITAL ROUTINE 'PEAK' 

PARAMETER (IFREQ=200) 

VIRTUAL lY (9000 ) 

DIMENSION INPUT( 15*IFREQ) 

DIMENSION OUTPUT( 10,70) 

DIMENSION ITABLE(68) ,VTYPE(2,2) 
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: :me:;s:om y( i450) 

COMMOM /?LOT/IX(2000) , Y 
-0YIV4LENCE ( IX, IMPUT) , (OUTPUT, Y( 501 ) ) 

DATA VTYEE/ ' VA ' , ' LLEY ' , ' BASE ' , ’ LIME ' / 

DATA ITABLE/1, 5, 5,800 ,1,0,0,61*0/ 

DATA IMLAST,IDIMO/600,70/ 

CALI £RRS£T( 112, .TRUE. , .FALSE. , .FALSE. , .FALSE. ,) 

D2REF=IDZDT0+( IDZDT1-IDZDT0)*.5 
C TYPE HERE IM VOL ’,NY,IY 

C SET TABLE VALUES FOR DIGITAL PEAK DETECTION ROUTINE 

C 

C ITABLE(6)=0 

ITABLE( 7 ) -0 
IMPTR=0 

T\tr 

MPlAKkS-0 

ITA5L£( 4)-100 !SET MIM DIFFERENCE FACTOR 

IF(IG.ATE.ME.O) THEN ! RESET IGATE AND IM IF REQUESTED BY USER 
I T .A 3 L £ ( 3 .) = I G AT ii 

y-iii 

END IF 
C 

c TRANSCRIBE DZDT DATA TO ARRAY INPUT 

C 

c N=0 

DO 1000 1=1, NY 
INPUT( I )=IY( 3000+1) 

C TYPE »,I,INPUT(I),IY(I) 

1000 CONTINUE 

C IADMIN=-16000 

IADMIM=IMIM( INPUT, NY ) 

C OPEN (UNIT=2,TYPE='NEW ,NAME=' INPUT.DAT' ,FORM='FORMATTED' ) 

C WRITE (2,1980) (IMPUT(I),I=1,NY) 

Cl 980 FORMAT(8I10) 

C CLOSE (UNIT=2) 

IF(IADMIN.LT.O) then ’OFFSET INPUT VALUES SO ALL ARE NON NEGATIVE 
DO 2000 1=1, NY 
INPUT(I)=INPUT(I)-IADMIN 
2000 CONTINUE 

END IF 

C TYPE *, (INPUT (I ) ,1=1, NY ) ,IADMIN 

C WRITE(2,*) (( INPUT (I )) ,1=1, NY ) , I ADMIN, SLOPE, INTRCP 

C 

C CALL ROUTINE TO FIND PEAKS AND TRAILING MINS TIMES 

C FOR dZ/dT MIN=PEAKS HEIGHT AND T=TIME AT TRAILING 

C MIMS - TIME WHEN dZ/dT=0 BEFORE PEAKS 

C OUTPUT ARRAY CONTAINS RESULTS OF PEAK DETECTION 

C 

CALL PEAK( ITABLE , INPUT, INLAST , INPTR, OUTPUT, IDIMO,NPEAKS) 

C TYPE *,NPEAKS, INLAST, INPTR 


50 


ORSGir^Ar !3 

OF POOR QUALITY 


c 


DED BECAUSE 


PEAK ROUTINE MAY AVERAGE INPUT FOR PEAKS 


DO 2500 1=1 ,MPEAKS 

OUTPUT { 2,1) = IN?UT( I FI X ( OUTPUT! 3.1))) 
CONTIi'iLiE 


IF(IADMIN.LT.O) THEN ! CONVERT RESULTS IF DATA OFFSET 
DO 3000 I=1,NPEAKS 
OUTPUT! 2 , I ) =OUTPUT! 2 , I )+I ADMIN 
OUTPUT! 4,1) =OUT?UT! 4 , I ) ADMIN 
OUTPUT! 7,1) =OUTPUT( 7 , I ) ADMIN 
30<I0 CONTINUE 

END IF 
C 

C REFINE PEAKS BASED ON SOME CRITERION IF NECESSARY 

C 

DO 4000 I=1,NPEAKS 

IF!OUT?UT! 2 , I ) . LE . DZREF ) GOTO 4000 

K?EAKS=KPEAKS +1 

IF!KPEAKS.EQ. I) GOTO 4000 

BO 3500 J=1 , 10 

OUTPUT! J , KPE A KS ) = OUT PUT ! J , I ) 

3500 CONTINUE 

4000 CONTINUE 

D WRITE (2, * ) NPEAKS, INLAST, INPTR 

D WRITE (2,20000) 

C WRITE O ,20000) 

20000 FORMAT!' PEAK NO . ' , 8X , ' AREA ' , 4X , ' P HEIGHT' ,6X, ' P TIME',4X, 

A 'L HEIGHT' ,6X, 'L TIME ',/, 1 1 X ,' HALF WIDTH', 4X,'T HEIGHT', 6X, 
B 'T TIME' ,8X, 'TYPE' ,8X, 'RATE'//) 

DO 4 L=l, NPEAKS 
KK=0UTPUT(9,L)+1 

D WRITE(2, 30000)!L,!0UTPUT!I,L),I=1,8),(VTYPE(M,KK),M=1,2), 

D 1 OUTPUT! 10, L) ) 

C WRITEd , 30000)!L,!0UTPUT!I,L),I=1,8).(VTYPE(M,KK),M=1,2), 

C 1 OUTPUT! 10, D) 

30000 FORMAT! 19, 5F 12. 0,/,9X,3F12.0,4X,2A4,F12.0) 

4 CONTINUE 

RETURN 
END 

FUNCTION IMIN ( INPUT, ICOUNT) 

DIMENSION INPUT! ICOUNT) 

IMIN=INPUT! 1 ) 

DO 1000 1=2, ICOUNT 
IF!INPUT!D .LT. IMIN) THEM 
IMIM=INPUT!I) 

END IF 

1000 CONTINUE 

RETURN 
END 

SUBROUTINE TFACT (NPEAKS, I FREQ, C, INLAST,K,TYMEAN, 

1 SLOPE, INTRCP,BSA,1Z0,IREC,IDZDT0,IADMIN,CATFIL,IDKREC) 
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PA.IAMETER (JEREQ-200) 

PARAMETER (MREC=15 ) 

PARAMETER (MUM=220) 

DIMEMSIOM lEPC(MUM) ,ZOEPC(MUM) 

DiMEMSIOM IMPUT( 15*JFREQ) .OUTPUT( 10 , 70) , DATCAR( 5 , 50) 

VIRTUAL IZO(9000 ) 

COMMON /PLOT/X( 1000) ,Y( 1450) 

REAL SL0PE(2) ,INTRCP(2) 

CHARACTER* 14 CATFIL 

EQUIVALENCE(Y(501 ) .OUTPUT) , (X, INPUT) , (DATCAR,Y( 1201 ) ) 

LOGICAL* 1 TYMEAM 
DATA I COUNT/ 0/ 

DATA lEPC/ 1,2, 3, 4, 

15,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24, 
125,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45, 
146,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66, 
167,63,69,70,71,72,73,74,75,76,77,78,79,80,81 ,82,83,84,85,86, 

1 87, 83, 89, 90, 91, 92, 93, 94, 95 ,96, 97, 98, 99,100,101,102, 

1103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118, 

M 19, 120, 121 . 122, 123, 12L, 125, 126, 127, 128, 129, 130, 131 , 132, 133, 134, 

2135.136,137,138,139,140,141,142,143,144,145,146,147,148,149,150, 

3151,152,153.154,155,156,157,158,159,160,161,162,163,164,165, 

4166.167,168,169,170,171,172,173,174,175,176,177,178,179,180, 

5181,182,183,184,185,186,187,188,189,190,191,192,193,194,195, 

6196,197,198,199,200,201,202,203,204,205,206,207,208,209,210, 

7211,212,213,214,215,216,217,218,219,220/ 

DATA ZOEPC/20.9,21 .0,21 .0,21 .0,20.9,21.0,21 .0,20.9,20.9,20.9,21.0, 

1 21 . 0 , 21 . 0 , 21 . 0 , 21 . 0 , 21 . 0 , 21 . 0 , 21 . 0 , 21 . 0 , 21 . 0 , 21 . 0 , 20 . 9 , 21 . 0 , 

1 21.0,20.9,20.9,21.0,21.0,20.9,20.9,20.9,20.9,20.9,20.9,20.9,20.9, 

2 20.9,21 .0,21 .0,21 . 1 ,22.4,22.4,21 .3,21 .1 ,21 .2,21 .2,21.2,21.2,21.2, 

3 21 . 3 , 21 . 2 , 21 . 2 , 21 . 2 , 21 . 2 , 21 . 5 , 21 . 3 , 2 1 . 2 , 21 . 3 , 21 . 4 , 21 . 3 , 21 . 2 , 21 . 3 , 

4 21. 3, 21. 3, 21. 3, 21. 4, 21. 3, 21. 5, 22. 2, 22. 2, 21. 4, 21. 1,21.3,21 .2,21.1, 

521.1,21.2,21.2,21.2,21.2,21.3,21.2,21.3,21.3,21.3,21.2,21.3,21.3, 
6 21.3,21.4,21.4,21.3,21.3,21.3,21.4,21.2,21.2,21.1,21.1,21.1,21.0, 
1 21 . 0 , 

721.1,21.0,21.2,21.3,21.3,21.3,21.3,21 .4,21.2,21.3,21.5,21.3,21 .4, 

8 21.3,21.4,21.5,21.4,21.4,21.5,21.5,21.3,21.2,21.2,21.1,21.2,21.1, 

9 21.1,21.1,21.1,21.1,21.4,21.4,21.4,21.3,21.5,21.3,21.3,21.3,21.5, 

1 21.4, 

1 21.4,21.4,21.5,21.3,21.3,21.3,21.3,21.7,21.3,21.1,21.1,21.7,21.9, 
1 21.3,21.3,21.1,21.8,21.1,21.2,21.3,21.1,21.0,21.2,21.2,21.2,21.2, 
1 21.2,21.2,21.2,21.3,21.3,21.3,21.3,21.1,21.2,21.2,21.3,21.3,21.1, 
1 21.1,21.3,21.2,21.7,21.7,21.2,21.2,21.2,21.1,21.2,21.2,21.2,21.2, 
1 21.5,21.3,21.5,21.2,21.3,21.3,21.2,21.2,21.3,21.1,21.2,21.2,21.2, 
1 21.2,21.4,21.2,21.2,21.2,21.2,21.3,21.3,21.3,21.4,21.2,21.31.3/ 

N=0 

C DO 200 I=1,NREC 

C KREC=IREC+I-1 

C PRINT *,IREC,KREC,N,NREC 

C CALL DIRECT ( KREC ,CATFIL , 2 , 1 , IZO) e- 

C DO 150 J=1,IFREQ OF POOR QUALITY 

C N=N+1 

C IZO (N)=ISPACE(J ,1) 

150 CONTINUE 
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OF POOR QUALITY 


100 

c 


IF( IADMIM.lt. 0) MDZDTO = IDZDTO-IADMIi‘ 
DO !00 1=1.2 

S L ' J ? I ) = ' 

iiTfriCP ( i ; = i . 

CONTINUE 
C=1 . 


4 00 


DO 1000 I=1,MPEAKS 
I3TRT=OUT?'JT(3.I) 

UOO •J = .STRT,2,-1 

TFdNPUT' J) .G-i’.NDZDTO ) GOTO 400 
TSTRT= J 
GOTO 500 
CONTINUE 
TSTRT=0 

UT'-r o^TjTT* to ’T.'^ >T\r >.Tr r^r;; T'cic 

IPTOLO FOR SFArtTIMG TIME 


LMc-uri. 


C L 2 5 C ON T I N ij E 

C IF( ICOUMT.GT. 1 ) THEM 

C DO 450 J=ISTRT,2,-1 

C IF( IFTOLD( J- I ) . LT . I FT0LD( J ) ) GOTO 4->0 

C TSTRT=600.*( ICOUNT-1 )-( 170. -J+1 ) 

C GOTO 500 

C450 CONTINUE 

D TYPE *, ’ ERROR MO BOTTOM ' 

C END IF 

C TSTRT=600.*(IC0UNT-1 ) 

500 CONTINUE 

TIME=(0UTPUT(8,I )-TSTRT) 

C CALL CGL(37,3,0) 

CALL CGL(33,TSTRT,350. ) 

C CALL CGL(37,5,0) 

CALL CGL(33,0UTPUT(8,I),350. ) 

ZMOT=IZO( I3TRT ) *SLOPE( 1 )+INTRCP( 1 ) 

C DO 550 J=1,MUM 

C IFdDKREC.EQ.IEPC(J)) THEN 

C ZNOT=ZOEPC(J) 

C GOTO 575 

C ENDIF 

550 CONTINUE 

575 CONTINUE 

DZDT=0UTPUT(2d)*SL0PE(2) + INTRCP(2) 

VNTVOL=C*TIME^( DZDT )/IFREQ/ZN0T**2 

DATCAR( 1 ,I)=VMTVOL 
IF(I.EQ.I) THEN 

BPM= IFRE0^60 . / ( OUTPUT (3.2) -OUTPUT( 3,1)) 

ELSE 

3?M= I FREQ^oO . / { OUTPUT (3d) -OUT.PUT( 3,1-1)) 

END IF 

C ELSE 

C END IF 
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c 


c 


c 

c 

c 

D 

D 

D 

D 

C 

c 


1 

1 

1 


600 

1000 


c 

c 

c 

C4000 


1000 


1500 


C 

C 

C 


CC-VNT70L*BPM/ 1 000 . 

CI-C0/3SA 

HI=0.0 

DATCAR(2, I)=C0 
DATCAR(3. I)=CI 
DATCAR(4,I)=DZDT 
DATCAR(5,I)=OUTPUT(3,I) 

SAVE TIME FACTORS FOR HI FROM DZDT DATA AND DZDT PEAK VALUES 
WHERE HI = DATCAR(4,I)/[Q-DATCAR(5,I)] 

AMD Q IS TIME OF Q WAVE BEFORE DZDT PEAK 


WRITE{2,* ) OUTPUT( 2 , I ) ,TIME, VMTV0L,TSTRT, J, ISTRT, 

ZMOT , C , DZDT , BPM , SLOPE , IMTRCP , IZO( ISTRT) 

WRITE (2,* ) (0UTPUT(3, I )-1 . )/IFREQ,VNTVOL,CO,CI ,HI , 

TIME/ IFREQ , DZDT , BPM , ZNOT , C 

PRINT I, TIME, ' ( ' ,0UTPUT(8,I), TSTRT, ' ) ',ZNOT, 

' ( ’ ,CUT?UT(3,I) , ' ) ',DZDT,' ( ’ ,0UTPUT(2 , 1 ) , ’ )' 
CONTINUE 
COMTIMUE 

OLDOUT =OUTPUT ( 3 , M PE AKS ) 

DO 4000 J=431 ,600 
IPT0LD(J-430)=IMPUT(J) 

CONTINUE 

RETURN 

END 

SUBROUTINE STRVOL ( KPEAKS , DATSD , IPEAKS) 

PARAMETER (IFREQ=200) 

DIMENSION DATCAR(5,50) ,DATSD(4),IPEAKS(KPEAKS) 

COMMON /PLOT/X( 1000) ,Y0 450) 

EQUIVALENCE (DATCAR , Y( 1201 ) ) 

DO 1500 1=1,4 

SUM=0.0 

SUMSQ=0.0 

DO 1000 J=1 , KPEAKS 
SUM=SUM+DATCAR( I , IPEAKS( J) ) 

SUMSQ=SUMSQ+DATCAR( I , IPEAKS( J ) ) **2 
CONTINUE 

DATCAR( 1,1) =SUM/KPEAKS 

DATSD(I)=SQRT(KPEAKS*SUMSQ-SUM*SUM)/KPEAKS 

CONTINUE 

RETURN 

END 

SUBROUTINE HICALC (NPEAKS, IFREQ, IECG,KREC,NREC,CATFIL) 
PARAMETER (JFREQ=200) 

VIRTUAL I ECO (9000 ) 

DIMENSION OUTPUT( 10,70) ,DATCAR(5, 50) ,INPUT( 15*JFREQ) 

CHARACTER* 14 CATFIL 

COMMON /PLOT/X( 1000) ,Y( 1450) 

EQUIVALENCE ( Y(50 i ) , OUTPUT) , ( Y( 1201 ) ,DATCAR) , (X , INPUT) 
N=0 


DO 200 1=1 ,NREC 
IREC=KREC+I-1 
TYPE *, CATFIL 


OF. POOR QUALITY 
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OF POOR QUALITY 


'3;’LL { [REC '^ATFIL,'* ^ TErn) 

EC- J-- 

J I 'J - . T 1 L IV_ y 

r ^ T _ ' T . 1 

I'CG (M) = [SPACE ( J ,l) 

■yO COi'-ITI;'.'L'E 

DO 1000 I=1,MPEAKS 
IEMD=DATCAR(5, I) 

i i: V, I . KE . 1 } THEM 

ISTRTr ( IEWD+DATCAR( 5 , I- 1 ) ) /2 

l-t r r' 

ISTRT:( l+IEMD}/2 
EMD IF 

ICOUMT=IEMD-ISTRT+1 

DO 500 J^1 , I COUNT 

I NPUT ( J ) = I ECG (I STRT+ J - 1 +6000 ) 

500 CONTINUE 

I T 'n'^ T M Y f T r _ 1 

T I ME = I TIME 

^ T t ^'~y <71 T I A A 

1 ( DATCAR (5,1) ~T IME ) / 1 FREQ 

IF(TIME) 600,600,700 
6O0» CONTINUE 

UhaCARvA,I)-0,0 
GOTO 1000 
700 CONTINUE 

DATCAR(4 , I ) ^DATCAR( 4 , I ) /TIME 
1000 CONTINUE 

RETURN 
END 

FUNCTION IMAX( INPUT, ICOUNT) 
DIMENSION INPUT( ICOUNT) 

IMAXr ICOUNT 

DO 1000 I=IC0UNT-1 , 1,-1 
IF(INPUT{I).GT.INPUT(IMAX)) THEM 
IMAX=I 
END IF 

1000 CONTINUE 

RETURN 
EMD 
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